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Executive  Summary 

Moeness  Amin  (PI) 


This  report  presents  the  results  of  the  research  work  performed  under  the  ONR 
funding,  grant  number  NOOOl 4-98- 1-0 176  over  the  period  of  November  15th,  1997  to 
September  30th,  1998.  The  research  efforts  on  this  grant  have  taken  four  different 
directions,  namely;  Array  Processing  Using  Time-Frequency  Analysis;  Fast  Computational 
Time-Frequency  Distribution  Kernels;  Blind  Spatial  Processing  for  Frequency  Diversity 
Spread  Spectmm  Communications:  Partial  Jamming  Suppression;  Adaptive  Array 
Processing  for  Multipath  Fading  Mitigation  Via  Exploitation  of  Filter  Banks.  In  the 
following,  we  highlight  the  approach  adopted  in  each  problem  and  summarize  the 
respective  results  and  contributions.  We  include  in  the  appendix  all  research  papers  which 
have  come  out  of  the  current  ONR  grant  since  its  inception  last  November.  Some  of  these 
papers  have  already  been  published,  while  others  are  either  under  review  or  will  be 
published  past  October  1st,  1998. 


1.  SPATIAL  AVERAGING  OF  TIME-FREQUENCY  DISTRIBUTIONS 

Time-frequency  distributions  have  been  recently  proposed  for  applications  to  array 
signal  processing  problems.  For  this  pu^ose,  spatial  time-frequency  distributions  (STFDs) 
have  been  introduced  and  represented  in  a  matrix  form.  The  elements  of  a  spatial  time- 
frequency  distribution  matrix  are  the  time-frequency  distributions  and  the  cross-time 
frequency  distributions  of  the  data  received  at  the  multi-sensor  array.  It  has  been  shown 
that  the  relationship  between  the  spatial  time-frequency  distributions  of  the  sensor  data  and 
the  time-frequency  distributions  of  the  sources  is  identical  to  that  of  the  sensor  data 
covariance  matrix  and  the  sources’  conrelation  matrix.  This  key  property  permits  direction 
finding  and  blind  source  separations  to  be  performed  using  tiie  sources’  time-frequency 
localization  properties. 

Blind  source  separation  consists  of  recovering  a  set  of  signals  of  which  only 
instantaneous  linear  mixtures  are  observed.  This  problem  has  been  typically  solved  using 
statistical  information  available  on  source  signals,  including  second  or  higher  order 
statistics.  Our  previous  contribution  to  this  area  has  shown  that  the  spatial  time-frequency 
distributions  is  an  effective  alternative  to  separating  sources  whose  signatures  are  different 
in  the  t-f  domain.  Successful  applications  of  STFDs  to  source  separation  require 
computing  STFD  at  different  time-frequency  (t-f)  points.  The  results  are  then  incorporated 
into  a  joint-diagnalization  technique  to  estimate  the  mixing,  or  the  array  manifold  matrix. 
This  matrix  transfers  the  source  signals  into  sensor  data  -  a  process  that  we  must  undo  by 
blind  separation,  often  up  to  a  multiplicative  complex  scalar  and  the  order  of  the  sources. 

We  have  continued  to  focus  our  research  work  on  array  processing  using  the  signal 
arrivals’  time-varying  frequency  characteristics.  We  have  concentrated  our  efforts  into  two 
different  specific  problems,  namely,  high  resolution  angle  of  arrival  estimation  and 
improved  blind  source  separations.  Both  problems  are  approached  using  the  signal  power 
distribution  over  two  joint  variables,  including  time  and  frequency. 

Our  first  contribution  is  improving  source  signal  estimation  by  performing  the  blind 
source  separation  using  ambiguity  functions  rather  than  time-frequency  distributions. 
Unlike  the  spatial  time-frequency  distributions,  where  the  separation  is  performed  using 
various  time-frequency  bins,  the  proposed  approach  is  based  on  information  from  the 
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ambiguity  domain.  The  procedure  is  similar  to  that  adopted  in  spatial  TFD-based  blind 
source  separation.  The  main  difference,  however,  lies  in  the  domain  in  which  the  mixing 
matrix  is  estimated.  In  the  proposed  technique,  the  mixing  matrix  is  obtained  by  joint 
diagonalization  of  the  spatial  ambiguity  (SAP)  matrices  rather  than  the  spatial  time- 
frequency  distribution  matrices.  These  SAF  matrices  are  formed  from  the  auto-  and  cross¬ 
ambiguity  functions  of  the  data  received  by  the  different  array  sensors.  It  is  shown  that  ^ 
each  time-lag  and  frequency-lag  point  in  the  ambiguity  domain,  the  sensors’  SAF  matrix  is 
related  to  the  sources’  ambiguity  function  matrix  in  the  same  fashion  as  the  sensors’  and 
sources’  TFDs.  Performing  blind  source  separation  in  the  ambiguity-domain,  rather  than 
the  t-f  domain,  provides  a  greater  ability  to  formulate  the  problem  using  the  signal  auto¬ 
terms.  Avoiding  cross-terms  both  simplifies  and  improves  the  performance  of  joint- 
diagonalization,  and  subsequently  enhances  source  estimation. 

Our  second  contribution  is  the  development  of  a  novel  approach  based  on  time- 
frequency  distributions  (TFDs)  for  separating  signals  received  by  a  multiple  antenna  array. 
The  sources,  as  stat^  above,  have  different  time-frequency  signatures  and  are 
instantaneously  mixed  at  the  array  sensors.  The  proposed  approach  provides  a  significant 
improvement  in  performance  over  the  recently  introduced  spatial  time-frequency 
distributions,  specifically  for  signals  with  close  time-frequency  signatures.  In  this 
approach,  spatial  averaging  of  the  time-frequency  distributions  of  the  sensor  data  is 
performed  to  eliminate  the  interactions  of  the  sources  signals  in  the  time-frequency  domain, 
and  as  such  restores  the  realness  property  and  the  diagonal  structure  of  the  source  TFDs, 
which  are  both  necessary  for  source  separation.  It  is  shown  that  the  proposed  approach 
yields  improved  performance  over  both  cases  of  no  spatial  averaging  and  averaging  using 
time-frequency  smoothing  kernels. 

Our  third  contribution  to  this  area  is  the  introduction  of  the  Time-frequency  MUSIC 
as  a  new  array  signal  processing  method  based  on  time-frequency  signal  representations. 
As  discussed  above,  spatial  time-frequency  distributions  have  successfully  been  used  to 
solve  the  problem  of  blind  source  separations  for  nonstationary  signals.  We  have  used  the 
same  underlying  structure  of  STFDs  to  solve  the  direction  finding  problem,  i.e.,  angles  of 
arrival  (AOA)  estimation.  A  new  method  for  the  estimation  of  the  signal  subspace  and  the 
noise  subspace  based  on  time-frequency  signal  representations  has  been  introduced.  The 
proposed  approach  consists  of  the  joint  block-diagonalization  (JBD)  of  a  set  of  spatial  time- 
frequency  distribution  matrices.  Once  the  signal  and  the  noise  subspaces  are  estimated,  any 
subspace  based  approach  can  be  applied  for  AOA  estimation.  We  propose  to  use  the 
MUSIC  algorithm.  The  effects  of  spreading  the  noise  power  while  localizing  the  source 
energy  in  the  time-frequency  domain  amounts  to  increasing  the  robustness  of  the 
eigenstructure  superresolution  method  with  respect  to  noise,  and  hence  improves  spatial 
resolution.  Performance  of  the  proposed  Time-Frequency  MUSIC  (TF-MUSIC)  for 
different  time-frequency  kernels  was  evaluated  numericily. 


2.  FAST  COMPUTATIONAL  TIME-FREQUENCY  DISTRIBUTION  KERNELS 

Time-frequency  distributions  (TFDs)  are  a  powerful  tool  for  nonstationaiy  signal 
analysis  in  a  wide  variety  of  applications.  However,  the  computational  burden  limits  the 
applicability  of  the  TFD  in  practice.  Hence,  design  of  fast  implementation  methods  are 
important  to  utilize  TFDs  for  applications  that  require  on-line  processing  such  as  radar  systems 
and  nonstationary  jammer  excision  in  spread  spectmm  communications.  The  TFD  is  often 
computed  using  the  following  procedure.  First,  tiie  bilinear  data  products  are  computed  which 
form  the  instantaneous  autocorrelation  function  (lAF).  Next,  the  lAF  is  smoothed  along  each 
time-lag  using  a  time-frequency  (t-f)  kernel  to  produce  the  local  autocorrelation  function 
(LAF).  Finally,  the  LAF  is  Fourier  transformed  along  the  time-lag  variable  to  yield  the 
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required  TFD.  Recognizing  the  fact  that  the  computations  required  to  compute  the  DFT  are 
much  smaller  than  those  needed  to  compute  the  LAP,  it  becomes  logical  to  target  the 
calculation  of  the  LAP,  for  significant  computational  savings. 

We  have  introduced  a  new  recursive  structure  for  efficient  computations  of  the  local 
autocorrelation  function.  This  recursion  is  achieved  by  decomposing  the  t-f  kernel  along 
each  time-lag  into  a  finite  number  of  cosinusoidal  terms.  Each  term  generates  one  recursive 
equation  to  update  the  LAP.  This  technique  is  named  the  trigonometric  decomposition  of 
time-frequency  kernels,  and  can  be  viewed  as  a  frequency  sampling  design  technique  of  t-f 
kernels.  This  proposed  class  of  Trigonometric  kernels  (also  referred  to  as  the  frequency 
sampling  (PS)  kernels),  while  providing  considerable  reaction  in  the  computation  of  the 
LAP,  retain  the  Hermitian  property  of  the  LAP  which  simplifies  the  DPT  operations. 
Purther,  trigonometric  decomposition  is  amenable  to  data  -dependent  kernel  design  and 
does  not  compromise  the  time-frequency  desirable  properties. 


3.  BLIND  SPATIAL  PROCESSING  POR  PREQUENCY  DIVERSITY  SPREAD 
SPECTRUM  COMMUNICATIONS:  PARTIAL  JAMMING  SUPPRESSION 

Prequency  diversity  spread  spectrum  (PD-SS)  has  been  recently  shown  to  be  a 
powerful  tool  for  digital  detection  as  well  as  an  effective  alternative  to  the  traditional  spread 
spectrum  techniques,  namely  direct  sequence  (DS-SS)  and  frequency  hopping  (PH-SS).  In 
a  general  context,  diversity  is  conceived  by  the  existence  of  several  replicas  (either  in  code, 
time  space,  or  frequency).  When  diversity  is  available  to  the  receiver  either  by  the  structure 
of  the  transmitted  signal  or  the  architecture  of  the  receiver,  optimum  spatial  signal 
processing,  which  is  blind  to  the  temporal  signal  characteristics,  can  be  derived. 

We  have  devised  a  novel  technique  to  obtain  optimum  blind  spatial  processing  for 
frequency  diversity  spread  spectrum  (FD-SS)  communication  systems.  The  sufficient 
statistics  for  a  linear  combiner,  which  prove  ineffective  due  to  the  interferer’s  spectral 
location,  are  modified  to  yield  improved  detection  under  partial  jamming  in  the  spectral 
domain.  Robustness  to  partial  time  jamming  is  achieved  by  extending  the  notion  of  replicas 
over  the  frequency  axis  to  repetition  over  the  time  variable.  Analysis  and  simulations  are 
provided  showing  the  advantages  of  using  FD-SS  with  spatial  diversity  to  combat 
interference  which  is  confined  to  a  narrow  frequency  or  time  support  relative  to  the  desired 
signal  extent  in  either  domain. 


4.  ADAPTIVE  ARRAY  PROCESSING  FOR  MULTIPATH  FADING  MITIGATION  VIA 
EXPLOITATION  OF  FILTER  BANKS 

The  spatial-temporal  equalization  can  be  achieved  by  space-time  adaptive 
processing  (STAP)  to  effectively  mitigate  inter-symbol  interference  (ISl)  and  co-channel 
interference  (CCI).  Such  a  scheme  consists  of  an  integrated  adaptive  array  and  a  temporal 
equalizer  to  perform  jointly  optimum  spatial  and  temporal  signal  processing.  However, 
solving  both  the  CCI  and  ISI  problems  simultaneously  by  conventional  STAP  methods  is 
difficult.  The  recent  STAP  methods  require  either  large  scale  matrix  inversion,  recursive 
calculation  or  cascaded  CCI  and  ISI  cancellers. 

We  have  proposed  an  efficient  subband  adaptive  array  processing  method  that 
utilizes  filter  banks  to  mitigate  both  the  CCI  and  ISI  effects  in  land  mobile  communications. 
A  subband  adaptive  array  has,  in  effect,  the  same  function  as  an  STAP  system,  while  the 
implementation  is  much  easier.  The  subband  adaptive  array  provides  sub-optimal 
performance  in  the  context  of  frequency-spatial  signal  processing,  which  enhances  the 
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signal  correlation  between  multipath  rays  prior  to  processing.  In  subband  adaptive  arrays, 
the  frequency  band  of  the  received  signal  is  divided  into  smaller  bands  through  the  use  of 
filter  banks.  Analysis  filters  yield  a  significant  increase  in  the  signal  coirelation  between 
the  multipath  rays  within  each  subband.  Such  an  increase  is  blind  in  the  sense  that  it  does 
not  require  a  priori  knowledge  of  the  arriving  signals.  As  a  result  of  increased  correlation, 
the  multipath  fading  associated  with  both  the  desired  and  the  interference  signals  is 
reduced,  thus  yielding  faster  convergence  of  the  adaptive  weight  vector  as  well  as  fewer 
degrees-of-freedom  (DOFs)  required  for  adequate  equalization.  The  importance  to 
properly  choose  the  analysis  and  synthesis  filters  to  obtain  good  equalization  performance 
by  subband  signal  processing  is  emphasized.  The  signal  correlation  enhancements  of 
multipath  signals  using  different  filter  banks  are  compared  and  several  simulation  examples 
are  provided  to  confirm  the  effectiveness  of  the  proposed  method. 


Research  Team  and  Collaborations 

The  research  team  working  on  this  grant  at  Villanova  University  are:  Professor 
Moeness  Amin  (Principle  Investigator),  Dr.  Yimin  Zhang  (Postdoctoral  Fellow),  Mr. 
Weifeng  Mu  (graduate  student),  and  Mr.  Govind  Mandapati  (graduate  student).  The  work 
by  the  two  graduate  students  have  not  matured  to  the  level  of  reporting  and  will  be  included 
in  the  next  year  report.  The  work  on  optimum  spatial  processing  and  partial  band  jamming 
is  performed  in  collaboration  with  Profi  Miguel  Lagunas  at  TSC  Department,  MODULO  5, 
Barcelona,  Spain.  The  work  on  fast  computation  time-frequency  kernels  is  joint  with  Dr. 
Gopal  Venkatesan  at  the  University  of  Minnesota.  The  time-frequency  MUSIC  is  the  result 
of  joint  collaboration  with  Dr.  Adel  Belouchrani  at  Ecole  National  Polytechnique,  Algeria. 
The  work  on  subband  arrays  is  in  conjunction  with  Dr.  Yimin  Zhang  while  he  was  in  Japan 
before  joining  Villanova  University  as  a  Postdoctoral  fellow,  in  April  1998,  and  his 
colleague  Dr.  K.  Yang. 
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SPATIAL  AVERAGING  OF  TIME-FREQUENCY  DISTRIBUTIONS 

Yimin  Zhang  and  Moeness  G.  Amin 

Department  of  Electrical  and  Computer  Engineering 
Villanova  University 
Villanova,  PA  19085,  USA 


ABSTRACT 

This  paper  presents  a  novel  approach  based  on  time-frequency 
distributions  (TFDs)  for  separating  signals  received  by  a 
multiple  antenna  array.  This  approach  provides  a  significant 
improvement  in  performance  over  the  recently  introduced 
spatial  time-frequency  distributions,  specifically  for  signals 
with  close  time-frequency  signatures.  In  this  approach,  spatial 
averaging  of  the  time-frequency  distributions  of  the  sensor  data 
is  performed  to  eliminate  the  interactions  of  the  sources  signals 
in  the  time-frequency  domain,  and  as  such  restore  the  realness 
property  and  the  diagonal  structure  of  the  source  TFDs.  which 
are  necessary  for  source  separation.  It  is  shown  that  the 
proposed  approach  yields  improved  performance  over  both 
cases  of  no  spatial  averaging  and  averaging  using  time- 
frequency  smoothing  kernels. 

1.  INTRODUCTION 

In  this  paper,  we  introduce  a  new  technique  for  source 
separation  based  on  time-frequency  distribution  methods.  The 
sources  have  different  time-frequency  signatures  and 
instantaneously  mixed  at  the  array  sensors.  The  number  of 
sensors  is  assumed  to  be  equal  to  or  greater  than  twice  the 
number  of  sources.  The  time-frequency  distributions  (TFDs)  of 
the  data  across  the  array  are  computed  and  used  to  construct 
spatial  time-frequency  distribution  matrices  (Situs).  By 
forcing  the  hermition  Toeplitz  structure  of  the  STFDs  and 
perform  spatial  symmetric  averaging  over  two  parts  of  the  array, 
we  achieve  significant  improvement  of  source  separation  over 
the  case  where  no  spatial  averaging  is  performed. 

Recently,  time-frequency  distributions  have  been  applied  to 
direction  finding  and  blind  source  separation  problems  in  array 
processing.  The  spatial  time-frequency  distributions  are 
introduced  in  [1]  and  represented  by  a  spatial  matrix  whose 
elements  are  the  time-frequency  distributions  of  the  data  across 
the  multi-sensor  array.  The  successful  application  of  STFDs  to 
separating  sources  with  identical  spectra,  but  different  tirne- 
frequency  signatures,  is  shown  in  [2].  In  this  application,  STFD 
matrices  computed  at  different  irf  points  are  incorporated  into  a 
joint-diagnalization  technique  based  on  generalized  Jacobi 
transform  to  estimate  the  mixing,  or  array  manifold,  matrix. 
This  matrix  is  then  used  to  estimate  the  sources’  signals  up  to  a 
multiplicative  complex  scalar  and  the  order  of  the  sources.  The 
general  theory  of  solving  blind  source  seperation  problems  using 
spatial  arbitrary  joint  variable  distributions,  including  those  of 
time  and  frequency,  is  given  in  [3].  In  [4],  the  two  arbitrary 
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variables  are  chosen  as  the  time-lag  and  frequency-lag,  and  the 
source  separation  was  performed  using  spatial  ambiguity 
functions.  The  use  of  STFDs  as  an  eigenstructure-based 
approach  for  direction  finding  is  given  in  [5],  where  the  Time- 
Frequency  MUSIC  technique  is  proposed  to  estimate  the  signal 
and  noise  subspaces. 

The  importance  of  joint-diagonalization  (JD)  in  the  STFD 
context  is  that  the  diagonal  structure,  the  distinct  eigenvalues, 
and  the  full  rank  properties  of  the  signal  TFD  matrix,  necessary 
for  source  separation,  can  be  easily  violated  when  operating 
with  a  single  t-f  point.  The  cross  time-frequency  distributions  of 
the  source  signals  yield  non-zero  complex  values  at  the  off- 
diagonal  elements,  rendering  the  estimation  of  the  mixing 
matrix  difficult,  or  even  impossible.  Also,  the  noise  contribution 
to  all  matrix  elements  at  low  SNR  cannot  be  ignored.  As  the 
interactions  of  the  source  signals  vary  over  the  time-frequency 
plane,  the  incorporation  of  several  STFD  matrices  at  different  t-f 
points  into  JD  enhances  diagonalization  and  leads  to  a 
successful  separation  of  signal  arrivals.  It  is  noted  that  the 
primary  motivation  of  using  smoothing  kernels  and  resorting  to 
other  variables  than  time  and  frequency,  specifically  the 
ambiguity-domain  variables,  is  to  allow  the  selection  of  joint- 
variable  points  where  the  interactions  of , the  source  signals  are 
insignificant. 

The  fundamental  role  of  the  proposed  technique  of  symmetric 
spatial  averaging  of  STFDs  is  the  effective  elimination  of  the 
signals’  intermodulations.  It  effectively  restores  the  diagonal 
structure  and  realness  property  of  the  signal  TFD  matrix. 
Symmetric  spatial  averaging  is  a  simple,  well-known  technique 
in  conventional  array  processing  [6].  It  uses  additional  array 
sensors  to  reduce  cross-correlation  in  coherent  and  correlated 
signal  environments,  and  thereby  permits  proper  angle-of-arrival 
estimations  and  source  separations.  It  is  shown  that  adopting 
this  technique  in  the  underlying  TFD-based  source  separation 
JD  problem  gives  robustness  to  t-f  point  selections  and  leads  to 
improved  performance  over  other  TFD-based  techniques, 
specifically  for  sources  whose  time-frequency  signatures  are  not 
very  distinct. 

2.  SPATIAL  TIME-FREQUENCY 
DISTRIBUTIONS 

The  data  vector  for  TV-element  array  is  given  by 

x(0  =  y(t)+n(t)  =  As(0+n(r).  (I) 

In  vector  forms,  x(r)=[xo(0»  •••»  is  a  noisy  instantaneous 

linear  mixture  of  the  source  signals  s{t)-[s\(t) . s„(t)Y  and 

n(r)  is  the  additive  noise.  The  mixing  matrix  A  is  the  transfer 
function  between  the  sources  and  the  array  sensors. 
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The  discrete-time  form  of  Cohen’s  class  of  TFD  for  signal 
x(t)  is  given  by  [7] 

DJtJ)  =  £  +  m  +  (2) 

where  t  and/ represent  the  time  index  and  the  frequency  index, 
respectively.  The  kernel  ^m,  /)  characterizes  the  TFD  and  is  a 
function  of  both  the  time  and  lag  variables.  The  cross-TFD  of 
two  signals  jc,  (r)  and  Xj(t)  is  defined  by 

('./)  =  £  OXt (t+m  +  l)x] (f + m - (3) 

/s-«o 

The  spatial  time-frequency  distribution  (STFD)  incorporates 
both  equations  (2)  and  (3),  and  is  defined  in  [2]  by, 

(r, / )  =  ^  (m,  /)  ®  x(r  -I-  m  +  /)x''  0  +  m  -  (4) 

/a— 

where  [!>«('./)]„=  (/,/),  for  i,  7=0,  N-1,  ®  denotes 

the  Hadamard  product,  and  is  the  time- 

frequency  kernel  associated  with  the  pair  of  the  sensor  data  JC/(r) 
and  Xj{t),  Under  the  linear  data  model  of  Eq.  (1)  and  assuming 
noise-free  environment,  the  STFD  matrix  takes  the  following 
simple  structure 

D„(r,/)  =  ADJf,/)A«  (5) 

where  T>ss{tj)  is  the  signal  TFD  matrix  whose  entries  are  the 
auto-  and  cross-TFDs  of  the  sources.  Eq.  (5)  is  similar  to  that 
commonly  used  in  conventional  blind  source  separation  and 
direction-of-arrival  (DOA)  estimation  problems  [8,9],  relating 
the  signal  correlation  matrix  to  the  data  spatial  correlation 
matrix.  If  T>ss{tJ)  is  a  full-rank  matrix,  the  two  subspaces 
spanned  by  the  principle  eigenvectors  of  T>xx{tf)  and  the 
columns  of  A  become  identical.  In  this  case,  directional  finding 
techniques  based  on  eigenstructures  can  be  applied.  If  is 

diagonal,  i.e.,  the  signal  cross-TFDs  at  the  time-frequency  point 
(f/)  are  zeros,  the  mixture  matrix  and  the  signal  waveform  can 
be  recovered  using  blind  source  separation  methods  [1,2].  In 
these  methods,  in  order  to  avoid  potential  problems  associated 
with  using  a  single  STFD,  STFDs  at  different  {tj)  points  are 
incorporated  into  a  joint-diagonalization  scheme.  Although  JD 
of  the  STFDs  is  effective  in  most  cases,  signals  with  close  time- 
frequency  signatures  are  still  difficult  to  separate.  As  shown 
below,  spatial  averaging  can  be  used  to  facilitate  signal 
separation. 

3.  SPATIAL  AVERAGING  TIME- 
FREQUENCY  DISTRIBUTIONS 

Symmetric  spatial  averaging  method  was  proposed  by  Pillai 
[6]  to  restore  the  full-rank  property  of  the  signal  covariance 
matrix  in  the  presence  of  coherent  signals.  In  this  section,  we 
extend  the  spatial  averaging  method  to  TFD  analysis,  and 
propose  the  signal  separation  method  by  joint  diagonalization 
(JD)  based  on  spatial  averaging  TFDs. 

Without  loss  of  generality,  we  consider  A/=2,  i.e.,  only  two 
sources,  ji(0  and  The  result  can  be  easily  extended  to 


multiple  sources.  By  ignoring  the  effect  of  noise,  the  received 
signal  at  /-th  array  sensor  is  represented  as 

x.{t)  =  xl^\t)-¥xl^\t)  =  5,  (6) 


where  G)i=27:sin(|)i/>.  (k  =1,2)  is  the  spatial  radian  frequency,  X  is 
the  RF  wavelength,  and  di  is  the  distance  between  0-th  and  i-th 
array  sensors.  We  assume  the  array  is  equi-spaced  linear  array. 
The  cross-TFD  of  x/(r)  and  Xj(t)  is 


J,  Sj  X,  Xj  X,  ij 

Jd  (tJ)  +  D 
L  J 

+\d  (,,/)  + D 

I  *2*2  *2 


X,  Xj 


itj) 


(7) 


Since  the  cross-terms  (second  term  in  each  bracket  in  (7))  are 
generally  complex,  it  is  clear  that  the  TFD  matrix  Dxx(//)  will 
not  provide  proper  phase  information  for  recovering  the  DOA  of 
the  arrived  signals  when  cross-terms  are  present.  However,  such 
phase  information  can  be  restored  by  using  spatial  averaging 
methods.  The  spatial  averaging  of  TFD  allows  the  signal 
separation  even  when  the  TFDs  of  multiple  signals  have  very 
similar  shapes  and  are  highly  overlapping. 

Let  the  number  of  array  sensors  be  2/V-lwith  the  array  center 
is  the  zeroth  sensor,  as  shown  in  Fig.l.  The  TFD  of  jco(0  and 
x,(f),i=0,  1,2,  ...,iV-l,is 

-K, 

where  we  note  ^o=0.  Similarly,  the  TFD  of  xo(0  and  X/(r)  is 

=  [d  (r,/)  +  D  (f./)  W'  +[z)„  +  D 

(9) 

The  spatial  averaging  of  (8)  and  (9)  is  given  by 


where 

fc,  =  D^(r./)+ReiD^,(r./)) 
fe,  =  D^,^(r./)+Re(D,^(r./)) 

Since  the  terms  in  the  brackets  are  all  real,  the  TFD  in  (10) 
correctly  represents  the  phase  information  caused  by  the 
propagation  delay  between  array  sensors,  even  when  the  cross¬ 
terms  are  complex.  The  matrix  formed  from  the  TFDs  (10) 
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is  hermitian  and  Toeplitz.  It  is  referred  to  as  the  spatial 
averaging  TFD  (SATFD)  matrix.  In  the  noise-free  environment, 
the  SATFD  matrix  can  be  expressed  as 

D„a,/)  =  AD.A''  (12) 

where 

D,(r,/)  =  diag[Z^,  b,]  (13) 

are  the  equivalent  TFD  of  the  signal  vectors.  Note  that 
no  longer  expresses  the  actual  TFD.  Clearly,  (12)  has  the  same 
format  as  (5),  and  D„(/,/)  here  is  diagonal  even  when  the 

cross-terms  of  the  TFD  of  the  signals  are  present.  Therefore,  the 
spatial  averaging  method  will  ensure  the  validity  of  the  TFD- 
based  signal  separation  in  the  presence  of  cross-TFD. 


4.  SIMULATION  RESULTS 


Equi-spaced  5-element  linear  array  is  used  for  simulation  with 
the  interelement  spacing  0.5X.  When  spatial  averaging  method 
is  used,  two  sub-arrays  are  formed,  each  with  3  elements.  Two 
sources  of  chirp  signals 


5,(0  =  5,(0  = 


(14) 


are  used,  where  \i  and  ©  are  chosen  to  be  O.OOStc  and  0.027:, 
respectively.  The  DOAs  of  the  two  signals  are  assumed  30°  and 
60°  from  the  broadside  direction.  No  noise  is  considered  here. 

Fig.2(a)  shows  the  Wigner-Ville  distribution  of  each  source 
signal,  and  Fig.2(b)  shows  the  respective  distributions  after 
signal  separation.  It  is  clear  that  the  array  fails  to  separate  5i(0 
and  52(0- 

In  the  TFD-based  signal  separation  method,  applied  in  Fig.  2, 
three  points  (t,f)  are  used  for  joint  diagonalization  at  /  =  32,  64, 
and  96.  The  frequency  /  is  chosen  so  that  the  TFD  at  the  first 
array  sensor  is  the  largest  for  a  given  t. 

To  show  the  effect  of  using  a  smoothing  kernel,  similar 
simulation  is  performed  with  the  Choi- Williams  kernel  [10] 
with  a=  0.1.  The  result  is  shown  in  Fig.3.  A  rectangular 
window  with  31  samples  in  both  time  and  frequency  scale  is 
used.  Since  the  two  signals  are  closely  spaced  in  the  t-f  domain, 
the  cross-terms  reduction  furnished  by  the  Choi-Williams  kernel 
is  limited,  and  again  the  array  fails  to  separate  the  two  signals. 

Fig.4  shows  the  separated  signals  under  the  same  conditions 
when  the  proposed  spatial  averaging  method  is  applied.  The 
signals  are  perfectly  separated,  except  for  their  order. 


5.  CONCLUSIONS 

Symmetric  averaging  of  spatial  time-frequency  distributions  has 
been  introduced.  The  averaging  improves  the  performance  of 
source  separation  using  joint-diagonalization  techniques.  It 
amounts  to  forming  a  spatial  hermition  Toeplitze  matrix  using 
the  time-frequency  distributions  of  the  data  across  one  half  of 
the  array.  This  matrix  is  then  added  to  the  spatial  matrix 
corresponding  to  the  other  half  of  the  array.  The  effect  of  this 
averaging  is  to  remove  interaction  between  the  source  signals  in 
the  time-frequency  domain.  Joint  digonalization  (ID)  using  a 


generalization  of  Jacobi  transform  is  then  applied  to  estimate  the 
mixing  matrix.  By  reducing  the  interaction  of  the  source  signals, 
the  JD  algorithm  yields  improved  performance  over  the  case 
when  no  averaging  is  performed.  The  paper  presented  an 
example  of  separating  two  chirps  signals  whose  time-frequency 
signatures  are  slightly  different.  The  proposed  approach  has 
successfully  separated  the  two  signatures,  while  other  non¬ 
averaging  methods  fail. 
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ABSTRACT 


In  previous  contributions,^”®  we  have  used  successfully  the  concept  of  spatial  time  frequency  distribution 
(STFD)  to  solve  the  problem  of  blind  source  separation  for  nonstationary  signals.  In  this  paper,  we  propose 
to  apply  this  same  concept  to  solve  the  problem  of  the  direction  of  arrival  (DOA)  estimation,  A  new  method 
for  the  estimation  of  the  signal  subspace  and  noise  subspace  based  on  time-frequency  signal  representation  is 
introduced.  The  proposed  approach  consists  of  the  joint  block-diagonalization  (JBD)  of  a  combined  set  of  spatial 
time  frequency  distribution  matrices.  Once  the  signal  and  the  noise  subspaces  are  estimated,  any  subspace  based 
approach  can  be  applied  for  DOA  estimation.  Herein,  we  propose  to  use  the  MUSIC  algorithm.  Performance 
comparison  of  the  proposed  approach  with  the  classical  MUSIC  method  is  provided.  The  influence  of  the  time- 
frequency  kernels  on  the  performance  of  the  newly  proposed  Time-Frequency  MUSIC  (TF-MUSIC)  is  evaluated 
numerically. 

Key  Words:  DOA  estimation,  Subspace  estimation,  Time-frequency  representation.  Spatial  time-frequency 
distributions,  Joint  block-diagonalization. 


1  Introduction 


In  many  signal  processing  applications,  signals  are  received  on  small  or  large  antennas  aperture.  In  some  of 
these  applications  the  multidimensional  signal  is  spatially  processed  and  reduced  to  one  dimension,  on  which  some 
desired  processing  are  applied.  This  often  occurs  when  dealing  for  example  with  time-frequency  signal  analysis.*^"® 
Typical  processing  in  this  case  is  the  signal  syntheses  from  the  time-frequency  plane.  In  other  applications,  the 
multidimensional  signal  is  directly  utilized  to  estimate  some  signal  parameters.  Subspace-based  methods  exploit 
this  spatial  information,  which  is  provided  by  the  multidimensional  signal.*^’®  These  methods  use  a  geometrical 
relation  involving  the  exact  moments  of  the  data,  commonly  the  data  covariance  matrix.  The  desired  signal 
parameters  are  then  extracted  by  solving  the  geometrical  relation  in  some  approximate  sense,  and  using  sample 
moments  instead  of  the  exact  ones.  This  class  of  methods  assumes  stationary  signals.  When  the  frequency  content 
of  the  measured  data  is  time-varying,  the  previous  approaches  can  still  be  applied.  However,  they  do  not  use  the 


*This  work  is  supported  by  ONR,  grant  #  N00014-98-1-0176 
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data  time*frequency  information,  which  if  properly  used  may  significantly  improve  performance. 

We  propose  to  develop  new  subspace-based  methods  with  a  geometrical  relation  involving  no  longer  the  exact 
moments  of  the  data  but  their  Spatial  Time  Frequency  Distributions.  The  spatial  time  frequency  distribu¬ 
tion  (STFD)  is  a  generalization  of  the  time  frequency  distribution  to  a  vector  signal.  In  previous  contributions, 
we  have  used  successfully  the  concept  of  STFD  to  solve  the  problem  of  blind  source  separation  for  nonstationary 
signals.^  In  this  paper  ,  we  study  the  application  of  this  new  concept  of  spatial  time  frequency  signal  representa¬ 
tions  to  solve  the  problem  of  DOA  estimation. 

This  paper  proposes  a  new  approach  based  on  a  joint  block-diagonalization  of  several  spatial  time  frequency 
distribution  matrices.  It  is  shown  that  performing  a  joint  block-diagonalization  of  a  combined  set  of  these  matrices 
provides  an  improved  estimate  of  the  signal  and  noise  subspaces  over  existing  techniques.  Once  the  signal  and  the 
noise  subspaces  are  obtained,  any  subspace  based  approach can  be  used  for  the  DOA  estimation.  Experimental 
performance  comparison  of  the  proposed  technique  with  the  classical  MUSIC^^  is  presented.  The  influence  of  the 
time-frequency  kernels  on  the  performance  of  the  proposed  Time-Frequency  MUSIC  (TF-MUSIC)  is  evaluated 
numerically. 

The  paper  is  organized  as  follows.  In  section  2,  the  data  model  is  stated  along  with  the  relevant  hypothesis. 
Spatial  time-frequency  distributions  are  recalled  in  section  3.  Section  4  presents  the  proposed  method  for  noise 
and  signal  subspace  estimation.  Section  5  discusses  the  general  conditions  that  should  be  satisfied  by  the  spatial 
time-frequency  distributions.  Numerical  simulations  evaluating  the  performance  of  the  new  method  are  given  in 
section  6. 


2  Data  Model  and  Assumptions 


Consider  an  array  of  m  sensors  receiving  waveforms  of  n  sources  (m  >  n).  The  data  snapshot,  representing 
the  sensor  output  vector  x(^),  is  assumed  to  obey  the  following  model: 

x{t)  =  A(^)s(t)  -f  n{t)  (1) 

where  n(f)  is  an  additive  noise,  6  =  [^i,  *  •  * , ^n]^,  A(0)  =  [a(^i ),•••, a(^„)]^,  and  s(i)  =  [si(^),  ■  •  • ,  is 

the  signal  vector  at  time  t,  with  the  superscript  ^  denoting  the  transpose  operator.  a(^/c)  is  the  transfer  vector 
between  Xk{t)  and  y(t).  The  parameter  vector  9  defines  the  DOAs. 

The  signal  vector  s{t)  is  assumed  to  be  a  nonstationary  multivariate  process  with 

1  ^ 

(i71)  :  R,  =  lim  —  ^s(f +  r)s(<)*  =diag[rii(r),---,r„„(r)]  (2) 

t=l 

where  superscript  denotes  the  conjugate  transpose  of  a  vector,  diag[-]  is  the  diagonal  matrix  formed  with  the 
elements  of  its  vector  valued  argument  and  rji(r)  =  limj.-j.oo  y  St=i  j  +  7-)s*(<)  denotes  the  autocorrelation 
of  Si{i).  Assumption  (HI)  means  that  the  component  Si(t),  1  <  «  <  are  mutually  uncorrelated  as  their 
cross-correlations  are  equal  to  zero. 

The  additive  noise  n(t)  is  modeled  as  a  stationary,  temporally  white,  zero-mean  complex  random  process 
independent  of  the  source  signals.  For  simplicity,  we  also  require  n(f)  to  be  spatially  white,  i.e. 

{H2)  :  E  {n{t  +  r)n*(t))  =  or^S{T)l  (3) 

where  <J(r)  is  the  Kronecker  delta  and  I  denotes  the  identity  matrix. 
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3  Spatial  time-frequency  distributions 


The  discrete-time  form  of  the  Cohen’s  class  of  time-frequency  distributions  (TFD),  for  signal  x(t),  is  given 

by^'  00  00 

DMtJ)=  E  E  <l>(m,l)x{t  +  m  +  l)x^it  +  m-l)e-^^^^‘  (4) 

1=^00  m=-oo 

where  t  and  /  represent  the  time  index  and  the  frequency  index,  respectively.  The  kernel  (!>{m,l)  characterizes 
the  distribution  and  is  a  function  of  both  the  time  and  lag  variables.  The  cross-TFD  of  two  signals  xi{t)  and 
X2{t)  is  defined  by 

E  E  +  m  +  +  (5) 

(=  —  00  m=  — CO 


Expressions  (4)  and  (5)  are  now  used  to  define  the  following  data  spatial  time-frequency  distribution  (STFD) 
matrix, 

CO  CO 

Dxx(t,/)=  E  E  4>(m,l)x(t  +  m-i-l)x*{t  +  m-l)e~^'^''^‘  (6) 

/=  — CO  m=  — CO 

where  [Dxx(^j/)]ii  ~  /)j  1)*’’?^* 

A  more  general  definition  of  the  spatial  time-frequency  distribution  matrix  can  be  given  as, 

CO  oo 

Dxx(t,/)=  E  E  (T) 

/=  — OO  m=— 00 

where  ©  designates  the  Hadamard  product,  and  [$(m, /)],_,•  =  <!>ij{m,l)  is  the  time-frequency  kernel  associated 
with  the  pair  of  the  sensor  data  Xi{t)  and  Xj{t). 

Under  the  linear  data  model  of  equation  (1)  and  assuming  noise-free  environment,  the  STFD  matrix  takes  the 
following  simple  structure; 

Dxx(t,/)  =  AWDss{<./)A(<  (8) 

where  Dss(U  /)  is  the  signal  TFD  matrix  whose  entries  are  the  auto-  and  cross-TFDs  of  the  sources  and  the 
superscript  ^  denotes  transpose  conjugation  for  a  matrix.  We  note  that  Dxx(t>  /)  is  of  dimension  mxm,  whereas 
Dss(t,/)  is  of  n  X  n  dimension.  For  narrowband  array  signal  processing  applications,  matrix  A{0)  holds  the 
spatial  information  and  maps  the  auto-  and  cross-TFDs  of  the  sources  into  auto-  and  cross-TFDs  of  the  data. 

Expression  (8)  is  similar  to  that  which  commonly  used  in  blind  source  separation^®  and  direction  of  arrival 
(DOA)  estimation  problems,  relating  the  signal  correlation  matrix  to  the  data  spatial  correlation  matrix.  The  two 
subspaces  spanned  by  the  principle  eigenvectors  of  Dxx(U/)  and  the  columns  of  A(0)  are,  therefore,  identical. 
Since  the  off-diagonal  elements  are  cross-terms  of  Dss(t,/),  then  this  matrix  is  diagonal  for  each  time-frequency 
(t-f)  point  which  corresponds  to  a  true  power  concentration,  i.e.  signal  auto-term.  In  the  sequel,  we  consider  the 
t-f  points  which  satisfy  this  property.  In  practice,  to  simplify  the  selection  of  auto-terms,  we  apply  a  smoothing 
kernel  (j>{m,l)  that  significantly  decreases  the  contribution  of  the  cross-terms  in  the  t-f  plane.  This  kernel  can  be 
a  member  of  the  reduced  interference  distribution  (RID)  introduced  in*"*  or  signal-dependent  which  matches  the 
underlying  signal  characteristics.^"’®'^® 
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4  Noise  and  Signal  Subspace  Estimation 


By  performing  the  SVD  of  the  steering  matrix, 

A(^)  =  [E,E„][D  of  V",  (9) 

and  incorporating  the  results  in  (8),  it  is  easily  shown  that 

Dxx{f.  /)  =  [E,E„]D(f,  /)[E,E„]",  (10) 

where  D(^,  /)  is  a  block-diagonal  matrix  given  by 

D(<,/)  =  diag[DV"Dss(<,/)VD  0],  (11) 

Since  E5  and  En,  which  respectively  span  the  signal  subspace  and  the  noise  subspace,  are  fixed  and  independent 
of  the  time-frequency  point  (t-f),  relation  (10)  reveals  that  any  matrix  Dxx(t,/)  is  block-diagonalized  by  the 
unitary  transform 

E=[E,  En].  (12) 

A  simple  way  to  estimate  the  signal  and  noise  subspaces  is  to  perform  the  SVD  on  a  single  matrix  Dxx(t,  /).  But 
indeterminacies  arise  in  the  case  where  Dss(t,  /)  is  singular.  To  avoid  this  problem,  we  propose  to  perform  a  Joint 
block- diagonalization  (JBD)  of  the  combined  set  of  {Dxx(</:5  /fc)|^  =  l../v  }  by  exploiting  the  joint  structure  (10) 
of  the  STFD  matrices.  This  joint  block-diagonalization  is  achieved  by  the  maximization  under  unitary  transform 
of  the  following  criterion, 

C'(U)  |u*Dxx(4,  /fc)u,P  (13) 

over  the  set  of  unitary  matrices  U  =  [ui,  •  — ,  Umj. 

The  BJD  criterion  allows  the  structural  information  contained  in  each  STFD  matrix  to  be  jointly  integrated 
in  a  single  unitary  matrix.  In  appendix  A,  we  propose  an  efficient  algorithm  for  solving  (13). 


DOA  estimation 

Once  the  signal  and  the  noise  subspaces  are  estimated,  one  can  use  any  subspace- based  technique  (e.g.  MU- 
SIC^°)  to  estimate  the  DOAs.  In  this  paper,  the  MUSIC  algorithm^®  is  applied  to  the  noise  subspace  matrix 
En  estimated  by  the  JBD.  We  obtain  the  Time- Frequency  MUSIC  (TF-MUSIC)  algorithm  which  estimates  the 
DOAs  by  finding  the  n  largest  peaks  of  the  localization  function, 

f{9)  =  \E^B.i0)\-\  (14) 


5  Conditions  on  the  time-frequency  kernels 

It  is  clear  from  (13)  that  from  a  perspective  of  the  joint  block-diagonalization  technique,  an  attractive  time- 
frequency  distributions  are  those  which  allow  the  formation  of  non-singular  diagonal  matrices  Dss(t,/)  at  the 
prospective  time-frequency  points.  This  property  should  be  viewed  in  light  of  the  following  observations: 


1.  The  off-diagonal  elements  of  Dss(<,  /)  should  be  zeros,  and  as  such,  distributions  that  mount  the  cross-terms 
on  the  top  of  auto-terms  should  be  carefully  studied  and  re-examined  for  qualification  under  the  proposed 
application  and  within  the  above  frame  work. 
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2.  Time-frequency  distributions  that  spread  cross-terms  over  the  entire  time-frequency  plane  should,  using  the 
same  argument,  lead  to  improved  performance  over  those  distributions  which  localize  the  cross- terms  in  the 
auto-term  regions. 

3.  Time-frequency  distributions  which  reduce,  but  still  localize  the  cross- terms  away  from  the  auto-term  regions 
appear  to  be  most  applicable  to  the  diagonal  matrix  requirements. 

4.  The  time-frequency  signatures  of  the  sources,  although  remain  distinct,  should  intersect  as  often  as  possible, 
producing  a  large  number  of  candidate  points. 

5.  Time-frequency  distributions  which  reduce  the  noise  variance  are  of  significant  importance. 

It  is  noteworthy  that  the  time-frequency  kernels  which  handle  cross-terms  differently  in  terms  of  their  smooth¬ 
ing  and  localizations  have  already  been  devised.  The  Wigner,^^  Choi- Williams, and  Born-Jordan  kernels^^  are 
good  examples  of  how  time-frequency  kernels  differ  in  cross-term  mitigations.  These  kernels  should  be  investi¬ 
gated  to  show  their  possible  distinct  offerings  in  the  proposed  method.  In  Section  6,  we  evaluate  the  performance 
of  the  new  technique  for  these  diiferent  kernels. 


6  Numerical  simulations 


Herein,  the  performance  of  the  proposed  TF-MUSIC  is  compared  to  the  MUSIC.^®  For  this  purpose,  we 
consider  a  uniform  linear  array  of  4  sensors  separated  by  half  a  wavelength  and  receiving  signals  from  2  sources. 
The  two  source  signals  arriving  at  0i  =  10^  and  60  =  -10^,  respectively,  are  unit  variance.  The  first  source 
signal  is  composed  of  a  chirp  signal  whose  frequencies  are  uji  =  O.Itt  and  cjo  =  O.Gtt.  The  second  source  signal 
is  composed  of  a  chirp  signal  whose  frequencies  are  =  O.Ctt  and  =  O.Itt.  The  noise  used  in  this  simulation 
is  zero-mean,  Gaussian  distributed,  and  temporally  white.  The  noise  power  cr-  is  adjusted  to  give  the  desired 
SNR=^  101ogio(^^“^)-  For  the  TF-MUSIC  algorithm,  50  STFD  matrices  are  considered.  The  variance  of  the 
estimated  DOAs  are  computed  over  100  independent  trials. 

Figures  1-4  illustrate  the  localization  of  the  two  sources  for  different  level  of  the  SNR  and  values  of  the  sample 
size.  In  this  simulation,  the  STFD  matrices  are  computed  using  Choi- Williams  kernels.  These  figures  shows  that 
TF-MUSIC  performs  better  than  MUSIC  for  low  SNR  level. 

Next,  we  compare  the  performance  of  TF-MUSIC  for  the  following  time-frequency  kernels: 


1.  Wigner  kernel 

2.  Choi- Williams  kernel 

3.  Born-Jordan  kernel 


and  with  respect  to  MUSIC. 

Figure  2  displays  the  variance  of  the  estimated  DOA  versus  SNR  for  500  samples.  The  solid  line  presents 
the  classical  MUSIC  algorithm.  The  dashed  line,  dashdot  line  and  the  dotted  line  correspond  to  the  TF-MUSIC 
using  Choi-Williams  kernel,  the  TF-MUSIC  using  Born-Jordan  kernel  and  the  TF-MUSIC  using  Wigner  kernel, 
respectively.  According  to  this  plot,  TF-MUSIC  using  Wigner  kernel  outperforms  the  other  methods.  These 
results  confirm  the  conditions  stated  in  Section  5  on  the  time-frequency  kernels. 
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Figure  1:  Source  localization  for  SNR=0  dB  and  100  samples. 


Figure  2:  Source  localization  for  SNR=10  dB  and  100  samples, 
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Figure  5:  Variance  of  6  vs  SNR. 


7  Conclusion 


In  this  paper,  a  new  DOA  estimation  approach  using  time-frequency  distributions  (STFD)  is  introduced.  The 
new  approach  is  based  on  the  joint  block-diagonalization  of  a  combined  set  of  spatial  time-frequency  distribution 
matrices.  The  later  are  made  up  of  the  auto-  and  cross-  TFDs  of  the  data  snapshots  across  the  multisensor 
array  and  they  are  expressed  in  terms  of  the  TFD  matrices  of  the  sources.  The  TFD  matrices  of  the  data 
and  sources  appear,  respectively,  in  place  of  the  spatial  and  signal  correlation  matrices  commonly  used  under 
stationary  environment.  The  diagonal  structure  of  the  TFD  matrix  of  the  sources  is  essential  for  the  proposed 
approach  and  is  enforced  by  incorporating  only  the  T-F  points  corresponding  to  the  signal  auto-terms.  The 
off-diagonal  elements  are  cross-terms  which  become  negligible  by  using  a  reduced  interference  distribution  kernel. 
We  have  focused  on  the  time-frequency  distributions  of  Cohen’s  class,  however,  one  can  use  any  other  bilinear 
time-frequency  signal  representations  such  as  the  affine  and  hyperbolic  classes. 

The  effect  of  spreading  the  noise  power  while  localizing  the  source  energy  in  the  time-frequency  domain 
amounts  to  increasing  the  robustness  of  the  proposed  approach  with  respect  to  noise.  The  paper  has  included 
numerical  experiments  of  simple  nonstationary  signals.  Performance  evaluation  of  the  proposed  method  for 
different  time-frequency  kernels  has  been  provided. 


Appendix  A: 

Joint  Block-Diagonalization  (JBD)  algorithm 


The  proposed  algorithm  for  the  joint  block  diagonalization  of  a  set  {M*  =  of  square 

matrices,  inspired  from,^^  maximizes  the  JBD  criterion  (13)  by  successive  Givens  rotations.  A  unitary  matrix 
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U  =  [uij]  is  sought  such  that 


Mj;.  =UM/,U^,  JSr  = 


(15) 


maximizes  criterion  (13). 


At  each  Givens  step,  we  choose  a  pivot  pair  p,q  and  make  Uij  =  Sij  except  for  the  elements  Upp  =  Uqq  =  cos(0), 
Upq  =  e^‘^sin(0)  and  Uqp  =  sin(0).  One  can  show  that  optimization  of  (13)  is  equivalent  to  the  maximization 

of  ^  ^ 

Q  *=  XI X!  I  sin(S)r7i^*^  |-  +  |  cos(6')m^p^  +  sin(0)mj*^  [-  (16) 

fe=i j=i 


for  p  =  1,  •  •  • ,  n  and  g  =  n  -1- 1,  •  •  • ,  m.  By  defining  the  vectors 


Q  also  reads 
where 


def 

[cos(^),sin(^)  cos(0),sin(0)  sin((^)] 

(17) 

glfe,- 

d^f 

K?. 

(k)  .  (kh 

(18) 

(k)  .  (Jfc)T 

(19) 

Q  =  v^Gv  =  3fie(G)  v  (20) 

G  =  X  X  +  g2fcjg2jj),  (21) 


and  we  have  used  the  fact  that,  G  being  hermit ian  by  construction,  its  imaginary  part  is  anti-symmetric  and  hence 
contributes  nothing  to  the  above  quadratic  form.  The  last  step  is  to  recognize  that  the  particular  parameterization 
of  V  is  equivalent  to  the  condition  v^v  —  L  Maximizing  a  quadratic  form  under  the  unit  norm  constraint  of  its 
argument  is  obtained  by  taking  v  to  be  the  eigenvector  of  3?e(G^G)  associated  to  the  largest  eigenvalue. 
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ABSTRACT 

Trigonometric  decomposition  of  the  time-frequency 
distribution  (TFD)  kernels  refers  to  representing  each  row  of 
the  kernel  in  the  time  time-lag  domain  by  a  finite  number  of 
cosinusoidal  terms.  This  allows  the  local  autocorrelation 
function  to  be  recursively  updated,  yielding  a 
computationally  efficient  TFD.  Unlike  the  spectrogram 
decomposition,  where  the  t-f  desirable  properties  are 
compromised,  the  proposed  trigonometric  decomposition 
preserves  both  the  support  and  marginal  properties  and 
permits  data-dependent  kernel  design. 

1.  INTRODUCTION 

Time-frequency  distributions  (TFDs)  are  a  powerful  tool  for 
nonstationary  signal  analysis  in  a  wide  variety  of 
applications  [1].  However,  the  computational  burden  limits 
the  applicability  of  the  TFD  in  practice.  Hence,  design  of 
fast  implementation  methods  are  important  to  utilize  TFDs 
for  applications  that  require  on-line  processing  such  as 
nonstationary  jammer  excision  in  spread  spectrum 
communications  [2]. 

The  discrete-time  Cohen’s  class  of  TFDs  is  expressed  as 
p(n,  l)x(m+ t)x'(m-  (1) 

fsT-i  m 

where  K  =  2L  +  1.  That  is,  the  TFD  can  be  computed  using 
the  following  procedure:  First,  the  bilinear  products 
x(m’h  t)x*(m~  i)2iTC  found.  This  two-dimenrional  function 
(function  of  time  and  time-lag)  is  called  the  instantaneous 
autocorrelation  function  (lAF).  Next,  the  lAF  is  smoothed 
along  each  time-lag  £  using  the  kernel  ^(rn£)  to  produce 
the  local  autocorrelation  function  (LAF).  Finally,  the  LAF 
is  Fourier  transformed  along  the  time-lag  variable  to  yield 
the  required  TFD.  The  Hermitian  property  of  the  LAF 
leading  to  the  realness '  of  the  corresponding  TFD,  is 
obvious  on  inspection  of  *(1). 

Early  attempts  at  devising  fast  implementation  of 
TFDs  used  the  Hermitian  property  of  the  LAF  to  provide 
efficient  frequency  transformations.  The  method  adopted  by 
Martin  and  FlandrinI3]  simply  computes  the  t-f  distribution 
as  the  real  part  of  the  complex  one-sided  FT  of  the  LAF. 
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Boashash  and  Black  [4],  on  the  other  hand,  have 
increased  the  throughput  by  combining  two  LAFs  in  one 
DFT,  and  doubling  the  number  of  output  frequency  samples. 
All  the  above  methods  are  exact,  as  they  introduce  no  error 
in  the  final  result.  But  they  serve  to  reduce  computations 
only  after  the  LAF  has  been  provided  with  cumbersome 
computations.  Hence,  they  are  potent  only  in  forming  the 
Pseudo  Wigner  distribution,  i.e.,  when  the  lAF  is  identical 
to  the  LAF. 

Spectrogram  decomposition  techniques  [5,6],  which 
directly  operate  on  the  data,  bypass  the  computation  of  the 
LAF.  They  compute  the  TFD  as  a  linear  combination  of  the 
spectrograms  (STFT)  of  the  data  with  different  window 
functions.  .  These  window  functions  are  the  dominant 
components  arising  out  of  a  spectral  decomposition  of  the 
kernel.  These  techniques  provide  approximations  of  the 
required  TFD  and  may  significantly  compromise  the  kernel’s 
desirable  time-frequency  (t-f)  characteristics  (e.g.  supports 
and  marginals)  unless  a  large  immber  of  decomposition 
terms  are  used. 

Recognizing  the  fact  that  the  computations  required  to 
compute  the  DFT  is  only  0(L  logL)  multiplications  whereas 
the  computation  of  the  LAF  is  O  (L^2]  multiplications,  it 
becomes  logical  to  target  the  calculation  of  the  LAF,  for 
significant  computational  savings.  A  new  class  of  recursive 
time-frequency(t-f)  kernels  that  can  be  used  to  recursively 
calculate  the  LAF  in  0(L)  multiplications  has  been  proposed 
in  [7].  This  approach  employs  modified  comb  filters  to 
average  the  bilinear  products  at  different  time-lags. 

In  this  paper,  we  generalize  the  recursive 
implementation  structure  presented  in  [7]  by  decomposing 
the  kernel  along  each  time-lag  into  a  finite  number  of 
.cosinusoidal  functions,  each  generates  one  recursive 
equation  to  update  the  LAF.  ‘  This  technique  is  named  the 
trigonometric  decomposition  of  time ’frequency  kernels^  and 
can  be  viewed  as  a  frequency  sampling  design  technique  to 
generate  computationally  efficient  kernels.  It  is  noteworthy 
that  this  proposed  class  of  kernels  (also  referred  to  as  the 
frequency  sampling  (FS)  kernels),  while  providing 
considerable  reduction  in  the  computation  of  the  LAF,  retain 
the  Hermitian  property  of  the  LAF.  Hence,  the  techniques 
outlined  in  [3-5],  can  be  used  in  conjunction  with  the 
recursive  computation  of  the  LAF  to  provide  additional 
computational  savings. 
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2.THE  FREQUENCY  SAMPLING  APPROACH 
We  first  briefly  review  the  frequency  sampling 
approach  for  filler  design  and  then  extend  it  to  the  TFD 
case.  The  FIR  frequency  sampling  design  method  uses  the 
DFT  to  produce  a  finite  impulse  response  filter.  This  design 
method  matches  the  desired  frequency  response  at  the  equally 
spaced  sampling  instants  and  approximates  the  desired 
response  between  the  sampling  instants.  Let 


*  =  1,2,..,// -1  (2) 
be  the  samples  of  the  desired  frequency  response  at  the 
sampling  instances  JtjQjj,  where  A/j^  and  Vj,  are  the  k-th 

sample  amplitude  and  phase  responses,  respectively  [8]. 
The  IDFT  of  H(k)  gives  the  corresponding  impulse 

response, ...  j 


/*  k=0 


(3) 


Setting  CCj^  =  ,  the  z-domain  transfer  function  of  the 


filter  can  be  expressed  as, 


H(k) 


.(4) 


1^-1  1  _  ^ 

=  S _ 

N  ^0  iv  jb.0 1  -  a*z 

It  is  clear  that  H(z)  has  a  recursive  form  even  though  the 
impulse  response  is  finite.  The  number  of  terms  is  equation 
(4)  is  equal  to  the  number  of  non-zero  frequency  samples. 
This  number  is  equal  to  the  number  of  frequency  samples  in 
the  passband,  providing  that  the  stopband  is  comprised  of 
zero  values.  Therefore,  filters  with  extended  stopband 
characteristics  lead  to  simplified  recursions.  It  is  this 
specific  property  that  makes  the  frequency  sampling  method 
computationally  attractive  for  TFDs. 

In  the  t-f  context,  the  frequency  sampling  filter  is  used 
to  process  the  data  bilinear  pTod’actsx(n-h  l)x*(n-‘  1)  and  is 
considered  for  each  time-lag.  Accordingly,  the  frequency 
samples  are  ambiguity-domain  variables  and  can  be  set  equal 
to  the  desired  values.  The  kernel  that  employs  these  filters 
at  different  lags  is  referred  to  as  the  frequency  sampling 
kernel.  In  the  case  that  all  desired  filters  have  zero-phase, 
the  kernel  becomes  real  and  symmetric,  and  is  given  in  the 
/  -  r  domain  by 

^  /  k^\  / 

p,=(N,-m 

By  taking  the  z-transform  over  the  time  variable,  then 

N,-i  Ni-ri 
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where  the  subscript /is  added  to  all  variables  to  indicate 
possible  dependence  of  the  parameters  on  the  lag  number. 


The  recursion  equation  to  update  the  LAF  r(n,l)  using 
the  t-f  kernel  in  (6)  is  given  by 

r(n.  1)  =  r{n  +  (N  - 1)/2. 1)  -  f(n  -(N  -  \)/2. 1).  (7-a) 

Where  F(nJ)  =  J^ri(n,l). 

k^O 


r^(n,  1)  -  Ci^r^in  - 1,  /)  -  r^(n  -  2, 1)  +  dt,[x(n  + 1) 
jc’(  n- l)-x(rt-l  +  l)x'(  n-l-l)], 

(7-b) 

/bfw,  1)  =  r^in  - 1,  /;  +  dQ,x(n  +  l)x’(n  - 1).  (7-c) 

and 

c„  =  2 cos(2kn/N,),  =  2(-lJ “M uCos(kK/ N ,)/ N 

^Ql  —  ^o/^l 

Time-frequency  kernels  using  expression  (5)  are  referred  to 
as  trigonometric  kernels.  A  large  stopband,  i.e.,  Mjy  =  0 

for  high  values  of  k,  make  them  computationally  attractive. 

To  satisfy  the  frequency  marginal,  we  should  have 
the  zeroth  frequency  sample  \j/(0,  l)=l.  Since  A/q/  “1* 
first  term  in  the  summation  (5-b)  requires  one  recursion 
multiplication.  For  the  time-support  property  to  be 
satisfied,  the  number  of  equally  spaced  frequency  samples 
at  the  time-lag  £  should  be  N ^  ^21  -hi*  For  a  cone- 

shaped  kernel,  Nf=2l‘hl*  This  condition  implies  that 

the  ambiguity-domain  frequency  'samples  must  be  more 
widely  spaced  for  smaller  lags.  As  such,  for  the  passband 
region  to  include  the  same  number  of  frequency  samples,  it 
must  be  narrower  for  higher  time-lags.  This  allows  recursion 
(7)  to  be  completed  at  the  same  time  for  all  lags  of  interest. 
We  maintain  that  this  specific  ambiguity  domain  passband 
characteristics  of  the  time-frequency  kernels  are  inherent  in 
all  low-pass  fixed  devised  kernels  including  the  Choi- 
Williams,  the  binomial,  and  the  Cone-shape  kernels  [1],  and 
is  satisfied  by  all  product  kernels  designed  for  reduced 
interference  distributions.  Fig  (1-a)  shows  the  two-term 
trigonometric  kernel  corresponding  to  using  two  frequency 
samples  in  the  passband.  Notice  the  pronounced  ripples  in 
the  passband  and  stopband.  This  situation  can  be  improved 
if,  instead  of  having  the  response  samples  go  directly  from 
unity  in  the  passband  to  zero  in  the  stopband,  we  allow  a 
transition  band.  Figure  (1-b)  shows  the  two-ierm 
trigonometric  kernel  obtained  when  the  second  passband 
sample  is  halved.  Clearly,  this  kernel  is  smoother  than  that 
shown  in  Fig.  (1-a). 

The  recursion  property  of  the  l-f  kernels  designed  b  y 
the  frequency  sampling  method  lays  the  ground  for  a 
computationally  attractive  kernel  decomposition  technique. 
By  applying  the  Fourier  series  expansion  to  any  given  t-f 
kernel  and  retaining  the  first  K  dominant  terms,  we  are,  in 
essence,  performing  a  frequency  sampling  kernel  design  in 
which  the  ideal  kernel  is  defined  by  the  K  frequency  samples 
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corresponding  to  those  terms.  In  turn,  existing  symmetric  t- 
f  kernels  can  be  made  amenable  to  recursive 
implementations  through  trigonometric  decomposition  in 
the  form  (3).  Therefore,  the  frequency  sampling  kernel 
design  method  and  the  trigonometric  kernel  decomposition 
are,  respectively,  the  design  and  the  analysis  sides  of  the 
same  problem,  and  each  can  be  used  to  formulate  recursive 
kernels.  It  is  noteworthy  that  the  one-term  trigonometric 
decomposition  of  all  kernels  satisfying  the  marginals  and 
the  time-support  properties' is  the  well  known  Bom-Jordan 
kernel  [1]. 

3.  SIMULATIONS 

To  illustrate  the  significance  of  the  trigonometric 
decomposition  of  t-f  kernels,  we  show  in  Fig.  (2-a)  the 
commonly  used  Choi-Williams  distribution  kernel  truncated 
to  satisfy  the  time-support  property.  The  two  term 
trigonometric  approximation  of  this  kernel  is  shown  in  Fig. 
(2-b).  It  is  evident  that  the  original  and  the  approximated 
kernels  are  very  much  alike,  but  the  trigonometric  kernel  is 
superior  in  computations. 

It  is  well  known  that  for  a  chirp  signal  the  autoterm  is 
not  located  near  the  a.xes  in  the  ambiguity  domain.  Hence, 
traditional  fixed  kernels  will  attenuate  the  chirp  autoterm  due 
to  their  inherent  lowpass  characteristics.  Here  we  consider 
the  specific  case  when  a  signal  consists  of  a  chirp  and  a 
sinusoid.  This  compounds  the  problem  further  as  it 
introduces  crossterms  that  lie  in  between  the  two  autoterm 
components  of  the  signal  in  the  ambiguity  domain.  A  data- 
dependent  kernel  designed  to  "match"  the  signal 
characteristics  will  have  values  close  to  unity  at  the 
autoterm  locations  and  values  close  to  zero  in  the  crossterm 
regions.  Such  a  data-dependent  kernel  can  be  simply 
constructed  by  placing  corresponding  unit/zero  values  at  the 
autoterm/crossterm  locations.  This  kernel  is,  in  general,  of 
infinite  extent  in  the  time-time  lag  domain.  From  the 
infinite  time  extent  kernel,  a  finite  time-extent  recursively 


implementable  kernel  can  be  constructed  using  the  FS 
method.  A  plot  of  this  designed  kernel  is  shown  in  Figure 
(3-a).  Note  that  this  is  the  minimum-variance  solution  for 
the  required  desired  signal  characteristics,  and  is  hence, 
statistically  optimal  for  white  noise [9]The  TFD  obtained  by 
the  designed  FS  kernel  is  shown  in  Figure  (3-b).  The  good 
performance  of  the  FS  kernel  in  depicting  the  t-f  signatures 
of  the  two  signals  is  evident. 
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ABSTRACT 

Blind  source  separation  consists  of  recovering  a  set  of 
signals  of  which  only  instantaneous  linear  mixtures  arc 
observed.  This  problem  has  been  typically  solved  using 
statistical  information  available  on  source  signals. 
Recently,  we  have  introduced  spatial  time-frequency  (t-f) 
distributions  as  a  new  and  effective  alternative  to 
separate  sources  whose  signatures  are  different  in  the  t-f 
domain.  This  paper  presents  a  new  blind  source 
separation  method,  exploiting  the  difference  in  the 
ambiguity-domain  signatures  of  the  sources.  The 
approach  is  based  on  the  digonalization  of  a  combined 
set  of  spatial  ambiguity  functions.  In  contrast  to 
existing  techniques,  the  proposed  approach  allows  the 
separation  of  Gaussian  sources  with  identical  spectral 
shape  but  with  different  ambiguity  domain  localization 
properties. 

1.  INTRODUCTION 

The  evaluation  of  quadratic  time  frequency 
distributions  of  the  data  snap  shots  across  the  array 
yields,  which  was  defined  in  [1]  as  spatial  time-frequency 
distributions,  permits  the  application  of  eigenstructure 
subspace  techniques  to  solving  a  large  class  of  channel 
estimation  and  equalization,  blind  source  separation,  and 
high  resolution  direction  of  arrival  estimation  problems. 
The  spatial  time-frequency  distribution  techniques  are 
most  appropriate  to  handle  sources  of  nonstationary 
waveforms  that  are  highly  localized  in  the  time- 
frequency  domain,  or  any  other  domain  of  a  different 
joint  variable  signal  representations. 

In  the  area  of  blind  source  separation,  the  spatial 
time-frequency  distributions  allow  the  separation  of 
Gaussian  sources  with  identical  spectral  shape,  but  with 
different  time  frequency  localization  properties,  i.e., 
different  signatures  in  the  time-frequency  domain.  On  the 
other  hand,  for  direction  of  arrival  estimation  problems, 
the  effects  of  spreading  the  noise  power  while  localizing 
the  source  energy  in  the  time-frequency  domain  amounts 
to  increasing  the  robustness  of  the  eigenstructure 
superresolution  methods,  which  are  based  on  both  signal 
and  noise  subspace  estimations,  with  respect  to  noise, 
and  hence  improve  resolution  performance. 
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In  this  paper,  we  present  a  new  approach  for 
blind  separation  of  sources  with  different  t-f  signatures. 
Unlike  the  spatial  time-frequency  distributions,  where 
the  separation  is  performed  using  various  lime-frequency 
bins,  the  proposed  approach  is  based  on  information 
from  the  ambiguity  domain.  The  procedure  is  similar  to 
that  adopted  in  spatial  TFD-based  blind  source 
separation.  The  main  difference,  however,  lies  in  the 
domain  in  which  the  mixing  matrix  is  estimated.  In  the 
proposed  technique,  the  mixing  matrix  is  abstained  by 
join  diagonalization  of  the  spatial  ambiguity  (SAF 
)matrices  rather  than  the  spatial  time-frequency 
distributions  matrices.  These  SAF  matrices  are  formed 
from  the  auto-  and  cross-  ambiguity  functions  of  the  data 
received  by  the  different  array  sensors.  It  is  shown  that 
the  sensor  spatial  ambiguity  function  matrix  is  related  to 
the  source  ambiguity  function  -matrix  in  the  same 
fashion  as  in  the  well  known  relationship  between  the 
sensor  and  source  covariance  matrices[2.3]. 

Performing  blind  source  separation  in  the 
ambiguity-domain  rather  the  time-frequency  domain 
provides  a  greater  ability  to  formulate  the  problem  using 
the  source  autoterms  rather  than  crosssterms.  Avoiding 
cross-term  points  in  any  joint  variable  representation 
leads  to  improved  source  separation,  as  it  will  have  a 
positive  effect  on  the  performance  of  the  joint- 
diagonalization  algorithm,  as  discussed  in  [1,4]. 

2.  TECHNICAL  APPROACH 

TFDs  are  a  powerful  tool  in  nonstationary  signal 
analysis[5.6].  So  far,  most  of  the  work  on  this  subject 
has  focused  on  temporal  signal  processing  without  much 
attention  given  to  the  spatial  variable.  In  the  following, 
we  introduce  a  new  approach  for  blind  signal  separation 
in  a  multisensor  array  using  the  ambiguity  function.  In 
order  to  present  the  spatial  ambiguity  function,  we  first 
recall  that  the  generalized  auto-  and  cross-ambiguity 
functions  of  two  signals Xjf/j  are  given  by 

D.je.'i)  =  x,(u+  r/2)x:(u-'r/2)e^^dH 

D,,(e.x)  =  ^y(e. i)J  x,(u+  T/2)x;(u-z/2)eJ^clii. 

(1) 
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for  ij=:l,2.  T  and  6  represent  the  time-lag  index  and 
the  frequency-lag  index,  respectively.  The  ambiguity- 
domain  kernel  Tj  is  a  function  of  the  above  two 
variables.  Expressions  (1)  are  now  used  to  define  the 
following  data  spatial  ambiguity  function  (SAP)  for  n 
sources  and  m  sensors 

D^(9, x)  =  (p(9,x) x(u  +  x/l)x*( It- x/2)e^^du{l) 
where 

lDJ9.T)lj  =  D,,  (e.T).  for  ij  =  \.2 . m. 

In  several  applications  such  as  semiconductor 
manufacturing  process,  factor  analysis,  narrowband  array 
processing,  and  image  reconstruction,  the  following 
linear  data  model  is  assumed, 

x(t)  =  As(t)  -h  n(t)  (3) 

where  the  m  x  n  matrix  A  is  referred  to  as  the  spatial 
matrix,  which  may  be  a  mixing  matrix  or  a  steering 
matrix,  depending  on  the  application  under 
consideration.  The  elements  of  the  m  x  1  vector  x(t), 
which  represent  the  measured  data,  are  weighted  sums  of 
the  source  signals,  and  as  such,  are  multicomponent 
signals.  On  the  other  hand,  the  elements  of  the  n  x  1 
vector  s(t)  are  source  signals  and,  therefore,  are  likely  to 
be  the  monocomponent  signals.  n(t)  is  an  additive 
noise,  often  assumed  zero  mean,  white  and  Gaussian 
distributed  process.  Due  to  the  linear  data  model, 

the  SAP  takes  the  following  structure 

Dje.x)  =  ADje.x)k‘<  (4) 

where  the  noise  is  neglected  as  a  first  step.  The 
superscript  “H”  denotes  the  complex  conjugate  transpose 
of  a  matrix  and  D^(6,x)  is  the  SAP  matrix  of  s(t).  We 
note  that  D^(9,r)  is  a  matrix  of  dimension  m-x-m, 
whereas  DJ9,t)  is  of  n-x-n  dimension..  Por 
narrowband  array  signal  processing  applications,  A 
holds  the  spatial  information,  and  through  a  similarity 
transformation,  it  maps  the  auto-  and  cross-APs  of  the 
monocomponent  signals  into  auto-  and  cross-APs  of  the 
data. 

Expression  (4)  is  similar  to  that  which  has  been 
commonly  used  in  blind  source  separation  and  direction 
of  arrival  (DOA)  estimation  problems,  relating  the 
signal  correlation  matrixUo  the  data  spatial  correlation 
matrix.  Here,  these  correlation  matrices  are  replaced  by 
spatial  ambiguity  function  matrices.  This  means  that  we 
can  solve  these  problems  in  various  applications  using 
the  new  formulation  which  is  more  tuned  to 
nonstationary  signal  environment. 

The  two  subspaces  spanned  by  the  principle 
eigenvectors  of  D^(9,x)  and  the  columns  of  A  are, 
therefore,  identical.  Since  the  off-  diagonal  elements  are 
cross-terms  of  D^J9,r),  then  this  matrix  is  diagonal 
for  allfftrjpoints  which  correspond  only  to  the  signal 
autoterms.  In  practice,  to  simplify  the  selection  of  such 


points  of  high  power  localization,  the  kernel 
that  significantly  decreases  the  contribution  of  the 
cross-terms  should  be  applied  and  the  ambiguity-domain 
points  near  the  vicinity  of  the  origin  should  be  selected. 

3.  JOINT  DIAGONALIZATION  USING  SAFs 
The  joint  diagonalization  for  blind  source  separation 
using  spatial  ambiguity  functions  consists  mainly  of  the 
following  two  steps: 

A.  Whitening 

Let  W  denotes  mxn  matrix  such  that 
f\VA)(VVAf  =  UU"  =  I.  i.e.,  WA  is  a  mxm  unitary 
matrix  (This  matrix  is  referred  to  as  a  whitening  matrix, 
since  it  whitens  the  signal  part  of  the  observations). 
Pre-  and  post-multiplying  the  spatial  ambiguity  function 
matrix  by  W  yields  the  whitened  SAP  matrix 

DJ0.r)  =  WDJ9.t)W"  (5) 

Prom  the  definition  of  and  equation  (5),  we  may 
express  T)  as 

DJ9.r)  =  UDJ9,x)U’'  (6) 

Since  matrix  U  is  unitary  and  D^(9,1:)  is  diagonal, 
expression  (6)  shows  that  any  whitened  data  SAP  matrix 
is  diagonal  in  the  basis  of  the  columns  of  the  matrix  U 
(the  eigenvalues  of  D^(9,X)  being  the  diagonal  entries 
of  D^(9,r)).  If,  for  the  (9^,r^)  point,  the  diagonal 
elements  of  D^(9^,rf)  are  all  distinct,  the  missing 
unitary  matrix  U  may  be  ‘uniquely*  (i.e.  up  to 
permutation  and  phase  shifts)  retrieved  by  computing  the 
eigendecomposition  of  However,  when  the 

0  -  T  signatures  of  the  different  signals  are  not  highly 
overlapping  or  frequently  intersecting,  which  is  likely 
to  be  the  case,  the  selected  (9^,T^)  point  often 
corresponds  to  a  single  signal  auto-term,  rendering 
matrix  D^(9f,Tf)  deficient.  That  is,  only  one  diagonal 
element  of  is  different  from  zero.  It  follows 

that  the  determination  of  the  matrix  U  from  the 
eigendecomposition  of  a  single  whitened  data  SAP 
matrix  is  no  longer  ‘unique’  in  the  sense  defined  above. 
The  situation  is  more  favorable  when  considering 
simultaneous  diagonalization  of  a  combined  set 

I /  =  !, . ,p}  of  p  SAP  matrices.  This 

amounts  to  incorporating  several  6  —  T  points  in  the 
source  separation  problem.  It  is  noteworthy  that  two 
source  signals  with  identical  0  —  T  signatures  can  not 
be  separated  even  with  the  inclusion  of  all  information 
in  the  9  -X  plane. 

B.  Joint  diagonalization  : 

The  joint  diagonalization(JD)  [7]  can  be  explained 
by  first  noting  that  the  problem  of  the  diagonalization 
of  a  single  n  x  n  normal  matrix  M  is  equivalent  to  the 
minimization  of  the  criterion  [8] 

C(M.V)=-J^  (7) 
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over  the  set  of  unitary  matrices  Hence» 

the  joint  diagonalization  of  a  set  \k  =  \,,k}  of  K 

arbitrary  n  x  n  matrices  is  defined  as  the  minimization 
of  the  following  JD  criterion: 

cm=-  ^C{M^  Vj  =  -X I V,  ‘m,v,  (8) 

k  ki 

under  the  same  unitary  constraint.  An  efficient  joint 
approximate  diagonalization  algorithm  exists  in  [7]  and 
it  is  a  generalization  of  the  Jacobi  technique  [8]  for  the 
exact  diagonalization  of  a  single  normal  matrix. 

Equations  (5-8)  constitute  the  blind  source 
separation  approach  based  on  SAP  which  is  summarized 
by  the  following  steps: 

•  Determine  the  whitening  matrix  Wfrom  the 
eigendecomposition  of  an  estimate  of  the  covariance 
matrix  of  the  data  (see  [7]  for  more  details) 

A 

•  Determine  the  unitary  matrix  U  by  minimizing  the 
joint  approximate  diagonalization  criterion  for  a  specific 
set  of 

whitened  SAP  matrices  j  I  /  =  1,..., p}t 

A  A 

•  Obtain  an  estimate  of  the  mixture  matrix  A  as  A 

A  0  A  ' 

=  W  17,  where  the  superscript  #  denotes  the  Pseudo¬ 
inverse,  and  an  estimate  of  the  source  signals  s(t)  as 

s(t)-^u  Wxity 

4.  SIMULATIONS 

In  Pigure  1,  we  show  an  example  of  the 
application  of  spatial  ambiguity  functions  to  the  blind 
source  separation  problem.  In  this  example,  we  have  two 
sources  and  three  sensors.  One  source  has  four  sinusoids, 
whereas  the  other  has  four  impulses.  The  TPDs  of  the 
two  waveforms  are  shown  in  Pig.  1-a.  The  TPDs  and 
APs  at  the  first  two  sensors  are  shown  in  Fig.  1-b. 
Using  the  spatial  ambiguity  functions,  we  are  able  to 
recover  the  original  signals  from  only  their  observed 
mixture,  as  shown  in  Fig.l-c.  It  is  clear  that  the 
recovered  signals  are  very  much  alike  the  transmitted 
ones,  which  shows  that  the  proposed  technique  undoes 
the  mixing  caused  by  the  propagation  channel. 


CONCLUSIONS 

A  new  blind  source  separation  technique  has  been 
introduced.  This  technique  is  based  on  the  ambiguity 
functions  of  the  data  received  by  the  different  sensors  of 
the  antenna  array.  The  proposed  technique  is  geared 
towards  the  applications  in  which  the  sources  have 
nonstationary  waveforms  that  are  localizable  in  the 
ambiguity  domain.  The  key  advantage  of  performing 
source  separation  using  spatial  ambiguity  functions 
rather  than  spatial  time-frequency  distributions,  as 
previously  suggested,  is  that  in  the  former,  the 
knowledge  of  the  signal  autoterms  being  positioned  near 
the  origin  simplifies  the  selection  of  the  appropriate 
signature  points  to  be  used  in  the  joint  diagonalization 
procedure  leading  to  signal  estimation. 
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TFDs  of  the  mixed  signals 


Ambiguity  functions  of  the  mixed  signals 

Fig.  1-b 
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ABSTRACT 

A  novel  technique  to  obtain  optimum  blind  spatial  processing  for  frequency  diversity  spread  spectrum 
(FDSS)  communication  systems  is  introduced.  The  sufficient  statistics  for  a  linear  combiner,  which  prove 
ineffective  due  to  the  interferer’s  frequency  characteristics,  are  modified  to  yield  improved  detection 
under  partial  jamming  in  the  spectral  domain.  Robustness  to  partial  time  jamming  is  achieved  by 
extending  the  notion  of  replicas  over  the  frequency  axis  to  a  repetition  over  the  time  variable.  Analysis 
and  simulations  are  provided  showing  the  advantages  of  using  FDSS  with  spatial  diversity  to  combat  the 
interference  when  it  is  confined  to  a  narrow  frequency  band  or  short  time  interval  relative  to  the  desired 
si^al  extent  in  either  domain. 
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1.  INTRODUCTION 


Frequency  diversity  spread  spectrum  (FDSS)  has  been  recently  shown  [1][2]  to  be  a  powerful  tool 
for  digital  detection  as  well  as  an  effective  alternative  to  the  traditional  spread  spectrum  techniques, 
namely  direct  sequence  (DSSS)  and  frequency  hopping  (FHSS)[3].  In  a  general  context,  diversity  is 
conceived  by  the  existence  of  several  replicas  (either  in  code,  time,  space,  or  frequency).  When  diversity 
is  available  to  the  receiver  either  by  the  structure  of  the  transmitted  signal  or  the  architecture  of  the 
receiver,  optimum  spatial  signal  processing,  which  is  blind  to  the  temporal  signal  characteristics,  can  be 
derived. 


In  this  paper  we  consider  the  spatial-frequency  diversity  problem  in  wireless  communication 
systems.  We  derive  optimum  spatial  processing  for  FDSS  systems  under  partial  band  jamming.  The 
principle  condition  for  the  optimality  in  these  systems  is  to  have  at  least  one  frequency  band  contaminated 
by  the  interference  and  another  band  that  is  jammer  free.  These  two  bands  will  here  in  be  refered  to  as  the 
‘hit’  and  ‘free’  bands,  respectively. 

In  partial  band  jamming,  the  above  condition  is  often  satisfyed,  as  the  desired  signal  is  present  in 
all  the  frequency  bands  in  which  the  transmission  bandwidth  is  slotted,  whereas  the  jammer  is  only  active 
in  few  bands  due  to  its  narrowband  frequency  characteristics.  In  this  case,  the  optimum  spatial  combiners 
can  be  obtained  in  a  two-step  design  procedure  as  follows.  First,  the  spatial  combiner  of  the  free  band  is 
computed  by  minimizing  the  mean  square  error  (MSB)  between  the  spatial  combiners’  outputs  of  the  free 
and  the  hit  bands.  In  order  to  avoid  trivial  solutions,  leading  to  null  beamformer  or  simultaneous  removal 
of  the  desired  signal  in  all  bands,  the  cross-correlation  of  the  two  band  outputs  is  set  equal  to  a  nonzero 
value.  Once  the  spatial  combiner  for  the  free  band  is  provided,  the  corresponding  output  signal  is  used  as 
a  time  reference  to  design  the  spatial  combiners  for  the  rest  of  the  bands.  It  is  shown  that  this  spatial 
design  procedure  is  blind  to  the  desired  signal  waveform  and  does  not  require  the  proper  labeling  of  the 
hit  and  the  free  bands  to  reach  the  optimum  solution. 
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The  sum  of  the  likelihood  of  the  different  bands  forms  the  sufficient  statistics  for  optimal  symbol 
by  symbol  detection.  This  requires  the  knowledge  of  the  jammer  level  in  every  hit  band  [1].  On  the  other 
hand,  the  sufficient  statistics  in  a  suboptimum  receiver  is  formed  only  the  bands  free  of  jammer  energy. 
The  performances  of  both  receivers,  however,  become  approximately  the  same  under  high  jammer-to- 
noise  ratio  in  the  hit  bands.  The  offering  of  spatial  processing  to  improved  detection  is  in  its  ability  to 
remove  the  jammer  contamination  in  the  hit  bands  prior  to  incorporating  them  into  the  detection  scheme. 
The  quality  of  jammer  suppression  in  the  hit  bands,  however  depends  on  the  angular  separation  of  the 
waveforms  received  by  the  different  antennas.  As  the  jammer  and  desired  signal  become  closely  spaced, 
removal  of  the  jammer  power  through  linear  combining  of  signal  arrivals  render  ineffective.  In  this  case, 
suboptimal  and  optimum  receivers  converge  to  the  same  solution.  It  is  noted  that  shadowing  occurs  when 
the  angles  of  arrival  of  the  jammer  and  the  desired  signal  coincide.  We  present  the  proper  mechanism  to 
form  a  sufficient  statistic  robust  to  the  shadowing  effect. 

Section  2  provides  the  general  stmcture  of  the  FDSS  systems  [5]  as  well  as  the  mathematical 
description  of  the  desired  signal.  Section  3  establishes  the  optimum  receiver  for  the  frequency  diversity 
spread  spectrum  signal  and  presents  the  decision  variable  and  probability  of  error  under  partial  band 
jamming.  In  section  4,  we  introduce  spatial  processing  for  improved  detection  in  FDSS  systems  and 
derive  the  optimum  solution  for  the  spatial  weights  in  the  hit  and  free  bands.  Section  5  describes  the 
optimum  receiver  robust  to  the  shadowing  effect.  Section  6  includes  an  adaptive  algorithm  as  an 
alternative  to  the  block  processing  described  in  section  4.  Section  7  reports  key  features  of  the  proposed 
design  procedure  for  optimum  spatial  processing.  It  is  shown  that  any  information  transmission  system 
using  coded  waveforms,  such  as  FDSS,  permits  blind  optimum  spatial  processing,  implying  that  the 
allocation  of  dedicated  time  or  frequency  slots  for  reference  framing  is  unnecessary.  Section  7  also 
includes  the  mechanism  for  exploiting  the  temporal  stmcture  of  the  desired  signal  to  blind  beamforming 
for  partial  time  jamming.  Finally,  a  set  of  simulations  illustrating  the  advantages  of  using  FDSS  combined 
with  spatial  diversity  are  provided  in  section  8. 
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2.  SPREAD  SPECTRUM  AND  FREQUENCY  DIVERSITY 


Direct  Sequence  and  Frequency  Hopping  are  the  most  familiar  techniques,  among  those  that  use 
spread  spectrum,  to  combat  bandlimited  Gaussian  interferers.  In  both  cases,  the  transmitted  signal  is 
controlled  by  properly  coding  the  carrier  waveform.  The  role  of  coding  is  to  enlarge  the  transmission 
bandwidth  either  by  fast  keying  in  the  time  domain  or/and  in  the  frequency  domain.  However,  the 
performance  of  FHSS  is  sensitive  to  partial  band  jamming  [2].  In  DSSS  systems,  high  levels  of  partial 
band  jamming  may  force  the  receiver  to  employ  a  whitening  filter,  together  with  a  sequence  detector  to 
remove  Inter-Symbol  Interference  (ISI)  caused  by  the  whitening  stage.  In  addition  to  the  highly  undesired 
ISI,  the  interference  excision  filters  increase  the  receiver  complexity  and  give  rise  to  self-noise,  which  is 
the  induced  correlation  across  the  PN  sequence.  To  mitigate  the  aforementioned  problems,  the  FDSS 
system  was  proposed  in  [1]  which  encompasses  the  principles  of  the  two  systems  DS  and  FH.  The  FDSS 
can  be  viewed  as  a  DS  scheme  where  successive  chip  frames,  within  a  symbol,  are  transmitted  in  parallel 
over  different  non  overlapping  frequency  bands.  On  the  other  hand,  the  major  difference  between 
traditional  FH  communication  systems  and  the  system  implementing  FDSS  is  that  in  the  later,  the 
information  symbol  is  repeated,  coded  or  uncoded,  in  several  frequency  bands,  instead  of  using  a  single 
band  at  time.  A  simple  mixture  of  both  systems  is  to  transmit  the  block  of  chip  symbols,  of  the  first 
scheme,  simultaneously  in  different  frequency  bands,  with  or  without  hopping.  The  flexibiliy  of  mixing 
the  DS  and  FH  schemes  grows  dramatically  when  a  single  information  symbol  produces  the  coded 
symbols  through  a  channel  encoder.  Since  all  possible  combinations  are  amenable  to  spatial  diversity 
processing,  the  transmitted  signal  architecture  depicted  in  Figure  1  will  be  used  to  introduce  the  proposed 
spatial  FDSS  scheme.  This  architecture,  which  was  used  by  several  authors  (see  [5]  as  an  example)  to 
encompass  all  SS  techniques  currently  reported,  allows  the  presentation  of  the  different  alternatives 
offered  by  channel  encoding,  waveform  coding,  and  frequency  diversity,  to  spread  the  spectrum  of  the 
transmitted  signal. 

The  first  stage  of  the  scheme  in  Figure  1  contains  the  channel  encoder.  Within  this  block,  the 
information  symbol  is  converted  into  coded  symbols.  The  coded  symbols  are  packed  in  blocks  of  Ni 


31 


symbols  and  these  blocks  are  arranged  in  N  branches.  Without  loss  of  generality,  it  is  assumed  that  all  the 
bands  have  the  same  bandwidth.  In  this  case  the  bandwidth  is  controlled  by  the  number  of  branches  N  as 
well  as  by  the  number  Ni  of  coded  symbols  per  branch.  This  presents  the  basic  contribution  of  the 
channel  encoder  to  the  bandwidth  structure  of  the  transmitted  signal.  When  the  number  of  symbols  per 
branch  is  equal  to  one,  and  the  coded  symbols  are  presented  in  the  N  channels  at  the  information  symbol 
rate,  there  will  still  be  a  coding  gain  at  the  receiver  even  when  using  a  pure  repetition  code.  This  type  of 
coding  amounts  to  the  lowest  complexity  case,  i.e.  Ni  is  equal  one  2md  the  information  symbol  is 
presented  simultaneously  to  all  branches  (pure  repetition  code  (N,l)). 

The  second  stage  is  the  waveform  coding.  Every  branch  is  multiplied  by  different,  or  equal,  chip 
signals  Yi  (t)  (i=l,N).  The  chip  rate  Xc,  as  a  fraction  of  the  coded  symbols  rate  Ts,  determines  the 
bandwidth  for  every  branch.  A  special  case  is  when  the  chip  rate  and  the  symbol  rate  coincide.  In  this 
case,  there  is  a  single  chip  symbol  Yni,i  per  branch,  and  no  spreading  is  produced.  Furthermore,  when  the 
alphabet  for  the  chip  symbols  is  (0,1),  the  transmitted  signal  can  be  viewed  as  associated  with  a 
Frequency  Hopping  system.  In  fact,  traditional  (single  band)  FH,  within  this  framework,  requires  a  single 
chip  is  to  be  on  per  a  symbol  interval ,  i.e.  a  single  branch  is  activated  at  a  time. 

The  third  stage,  the  last  in  this  sequence,  contains  the  pulse  shaping.  The  shaping  functions  \|/i(t) 
have  different  frequency  support.  Multicarrier  CDMA  or  OFDM  can  be  encompassed  in  the  proposed 
scheme,  depending  on  the  design  of  this  last  stage  in  the  overall  structure  of  signal  generation.  Assuming 
non-overlapping  frequency  support  for  the  shaping  functions  produces  a  transmitted  signal  which  is 
usually  referred  to  as  Frequency  Diversity  spread  spectrum  signal.  The  non-overlapping  frequency  bands 
often  have  the  same  bandwidth  Bq.  This  bandwidth,  which  includes  the  roll-off  factor,  has  to  be  the 
minimum  bandwidth  expected  from  the  jammer.  This  choice  guarantees  that  no  narrowband  intereferer  is 
present  for  every  frequency  diversity  band;  otherwise,  a  whitening  filtering,  together  with  a  Viterbi 
equalizer  to  remove  the  ISI,  will  be  needed.  The  filter  bank  is  implemented  as  a  polyphase  filter  [6]. 
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Finally,  note  that  within  a  FDSS  band,  the  system  may  have  independent  frequency  hopping  [4], 
allowing  phase  continuity  and  coherent  optimal  combining  of  the  received  replicas. 

The  formulation  of  the  transmitted  signal  is  shown  in  (1)  for  the  general  case  when  code,  time, 
and  frequency  diversity  are  used  simultaneously  to  spread  the  spectrum.  In  (1)  n(.)  denotes  the 

rectangular  pulse,  Nc  the  number  of  chips  per  coded  symbol  interval,  Tj  the  symbol  duration,  Nj  the 
number  of  coded  symbols  in  every  branch  per  information  symbol,  and  Xc  the  chip  duration. 


N  Nc 

XT(t)=X  Z  Z  X  ai.n.m  11 

im=\n=l  /=1 


^t-iTs -niTs/NQ-Itc  ^ 


V 


TsKNi*Nc) 

In  the  case  of  no  spreading,  Ni=l,  and  equation  (1)  simplifies  to  (2). 


Yi,n.m,l  Vm(t"iTs) 
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N  Nc 
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Finally,  for  a  single  chip  per  symbol  interval  and  the  repetition  code,  the  transmitted  signal  is  (3). 


(2) 


Xrlt)  = 

m  l=:l 

This  signal,  referred  to  as  FDSS  signal,  will  be  used  from  now  on  without  loss  of  generality. 


(3) 


3.  OPTIMUM  RECEIVER  FOR  FDSS 

This  section  highlights  the  differences  between  optimum  and  suboptimum  detection  in  terms  of 
performance  versus  complexity.  It  includes  results  from  [1]  which  are  most  pertinent  to  the  underlying 
problem  and  which  serve  the  derivation  of  the  optimum  receiver,  including  spatial  diversity,  in  the 
following  sections. 


The  decision  variable  Ai  is  a  sufficient  statistic  to  perform  optimum  symbol  by  symbol  detection. 
This  decision  variable  contains  the  following  terms:  the  received  samples,  after  matched  filtering  and 
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sampling  at  the  symbol  rate,  Zi,m ;  the  chip  symbol,  assumed  to  be  known  by  the  receiver  Yi,m ;  and,  factors 
Fm  which  depends  on  the  spectral  density  of  the  noise  plus  the  jammer  in  every  frequency  band  m. 


N 

m=l 


(4) 


Since  the  bandwidth  of  every  band  is  equal  or  less  than  the  jammer’s  bandwidth,  the  jammer  and  the 
noise  are  assumed  to  have  flat  spectral  densities,  given  by  Jo  and  No  respectively.  In  this  case,  the  factor 
Fm  in  the  above  equation  is  given  by  (5). 


F„=l/(1+ 

Fm  =  1  for  free  bands 

The  decision  for  the  received  symbol  is  the  alphabet  symbol  a;  which  maximizes  Re  (  ai .  Ai ), 
where  Re(.)  indicates  the  real  part  of  a  complex  argument.  The  resulting  probability  of  a  symbol  error  is. 


A 

N„ 


for  bands  hit  by  the  jammer 


(5) 


Pe=Q 


2^ 

Nr, 


(6) 


where  Ej  is  the  received  energy  per  symbol,  Q(.)  is  the  error  function,  and  p  is  the  loss  factor  which 
reflects  the  degradation  due  to  the  presence  of  jammer.  This  factor  depends  only  on  the  spectral  density  of 
the  jammer,  normalized  by  the  noise  density,  as  well  as  the  number  of  hit  and  free  bands,  denoted  by  Nh 
and  Nf  respectively,  relative  to  the  total  number  of  frequency  bands  N. 


Nf 


(7) 


The  above  factor  can  be  modified  accordingly  in  the  case  of  unequal  jammer  strength  for  every 
band.  When  the  jammer  has  the  same  strength  in  all  the  hit  bands,  the  above  formula  can  be  expressed  as 
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a  function  of  the  fraction  of  bands  hit  H-  This  factor  is  the  key  parameter  in  the  performance  of  FDSS 
systems  [2]. 


(i+r« 


O-j?) 


A  suboptimum  receiver  does  not  use  the  side  information  corresponding  to  the  jammer  strength, 
and  the  decision  variable  A  is  only  formed  by  the  free  bands.  The  probability  of  symbol  error  in  such  a 
case  is  given  by  (9),  which  shows  that  the  loss  factor  P  is  equal  to  1-T|,  being  independent  of  the  jammer 


strength. 


Note  that  both  receivers  degrade  in  performance  when  gets  close  to  one;  in  other  words,  when 
the  jammer  spreads  its  bandwidth  covering  the  entire  transmission  band.  A  detailed  study  of  both 
receivers  and  their  comparison  with  DS  and  FH  systems  can  be  found  in  [1,2,4],  showing  the  advantages 
in  many  respects  of  FDSS  over  traditional  DSSS  and  coherent,  or  incoherent,  FHSS  for  partial  band 
jamming.  It  should  be  noted  that  for  a  DSSS,  the  processing  gain  is,  in  essence  l/q,  and  the  equivalent  p 

is  given  by: 


which  is  always  lower  than  that  of  in  the  FDSS  optimum  receiver. 

Spatial  diversity  in  FDSS  is  next  introduced  in  order  to  alleviate,  whenever  it  is  possible,  the 
impact  of  the  hit  bands  on  the  receiver  performance.  This  is  to  say  that  spatial  processing  strives  to  reduce 
the  effective  T|,  by  removing  jammer  contamination  in  the  hit  bands. 
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4.  SPATIAL  PROCESSING  IN  FDSS  SYSTEMS 


As  mentioned  in  the  previous  Section,  when  the  number  of  bands  hit  by  the  interference  is  close 
to  N,  the  total  number  of  frequency  bands,  both  the  optimum  and  suboptimum  receivers  degrade  in 
performance.  Furthermore,  high  levels  of  jamming  tend  to  impair  the  performance  of  both  receivers, 
rendering  all  efforts  of  labeling  every  band  useless.  As  it  is  typical  in  digital  transmission,  when  temporal 
processing  is  insufficient  to  counteract  interferences,  spatial  diversity  may  be  introduced  to  remove,  or  or 
at  least  reduce,  the  jammer  contamination  of  the  desired  signal,  independent  of  its  temporal 
characteristics.  Taking  advantage  of  the  frequency  diversity  available  in  FDSS  systems,  we  propose  a 
procedure  to  design  optimum  spatial  diversity  processing  that  is  blind  to  the  desired  signal  waveform.  In 
fact,  the  procedure  is  valid  for  all  the  SS  schemes  which  implement  the  architecture  of  Figure  1  for  the 
transmitted  signal  generation,  providing  N  is  greater  than  one. 

It  is  assumed  that  the  desired  signal  is  present  in  all  the  frequency  bands.  It  is  also  assumed  that 
the  labeling  of  hit  bands  is  properly  done  (This  assumption  is  latter  relaxed).  All  the  air  interfaces,  down 
conversion,  and  filtering  channels  are  calibrated.  Full  coherence  is  assumed  for  the  desired  signal  across 
spatial  diversity  channels. 

Let  us  assume  that  2£f.nand  2Li,nare  the  snapshots,  at  time  n,  of  the  free  and  hit  bands,  respectively. 
Both  bands  are  selected  to  start  up  the  design  of  the  respective  optimum  spatial  combiners  ^  and 

Once  the  chip  symbol  is  removed  for  the  two  bands  selected,  i.e.  2£f.n  =>  7/,„-2£f.n  and  2^,„  „ .  2^,„ . 

the  desired  signal,  which  is  present  in  both  bands  can  act  as  the  reference  signal.  The  design  of  both 
combiners  can  then  be  performed  through  the  quadratic  cost  function  minimization: 


^  =  E 


(11) 
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The  undesired  solutions  to  the  above  objective  are  either  the  null  vector  for  both  combiners,  or  the  case 
when  the  desired,  together  with  the  jammer,  is  nulled  out  simultaneously.  A  suitable  technique  to  avoid 
both  solutions  is  to  set  the  cross-correlation  of  both  combiner’s  outputs  to  some  chosen  value. 


where  ^  is  some  constant  different  from  zero 

Naming  the  array  autocovariances  and  cross-covariance  of  the  two  frequency  channels  as, 

(13) 


the  minimization  problem  becomes: 

Minimize  ^  =  }Vf  +w^  R^Wf^-w^  Pf  ^Wf  (14.a) 

Constrained  to  ^  — / /?— /  ”  ^  (14.b) 

Forming  the  Lagrangian  and  setting  partial  derivatives  equal  to  zero,  the  optimum  spatial  combiners  are 
obtained  as: 

IS*  ('5'>) 

where  X  is  the  Lagrange  multiplier.  Since  %  is  minimum  for  minimum  X,  the  optimum  combiner  for  the 
jammer  free  band  is  the  eigenvector  of  (15.a)  associated  with  the  minimum  eigenvalue.  After  v^is  found, 
the  combiner  for  the  hit  band  is  just  the  Wiener  solution  for  the  cross-correlation  vector  Once 
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the  combiner  of  the  free  band  is  derived,  the  rest  of  the  combiners  can  be  derived  by  using  the  band  free 
output  as  a  time  reference  for  all  other  bands.  The  unconstrained  minimization  procedure  is  shown  in 
(16.a),  together  with  the  solution  for  the  combiners  in  (16.b).  Since  every  frequency  band  has  to  be 
calibrated  in  the  front-end,  the  reference  combiner  is  normalized  for  response  equal  to  one  in  a 
preselected  diversity  channel  (i.e.  Wf(l)=l)  before  it  is  used  in  (16.b);  at  the  same  time  this  normalization 
also  helps  the  weights  implementation  with  finite  register  length. 


Xn.f  -»:? 


H 


min.  for  /2  9^  / 


-lOT/n 


Wu  =  A  ,  r  ,  -w  f 

-h  =h,f-J 


(16.a) 


(16.b) 


It  is  noteworthy  that,  since  the  jammer  is  uncorrelated  from  one  band  to  another  due  to  the  chip 
demodulation,  any  two  bands  can  be  selected  to  start  up  the  design  in  (16.a).  Once  a  combiner  is  derived, 
its  output  provides  a  time  reference  for  the  rest  of  the  bands.  This  implies  that  labeling  is  not  essential  to 
design  the  spatial  processor.  Nevertheless,  when  required,  a  suitable  procedure  to  determine  the  free  and 
the  hit  bands,  to  be  used  in  (15. a),  is  to  choose  the  two  bands  with  minimum  and  maximum  received 
power  respectively. 


Note  that  when  the  transmission  bandwidth  is  very  small  compared  with  the  central  frequency, 
and  adequate  labeling  of  free  and  hit  bands  is  available,  there  is  a  great  saving  in  computations  since  the 
combiners  derived  in  (15),  for  a  pair  of  free  and  hit  bands,  can  be  applied  directly  in  the  rest  of  the  bands. 


To  gain  insight  into  the  design  method,  we  inspect  the  design  equation  for  the  case  when  the 
desired  signal  is  a  point  source;  in  this  case,  the  cross-covariance  matrix  is  rank  one. 


the  corresponding  equations  for  the  optimal  combiners  reduce  to 


(17) 
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as  expected  from  the  optimum  combiner  equations,  when  the  vectors  of  the  second  term  of  (17)  are 
known.  Nevertheless,  the  rank  one  approximation  of  the  cross-covanance  provides  poor  performance 
when  the  number  of  snapshots  used  in  its  estimation  is  small.  In  general,  the  formulation  (15)  offers  much 
better  results  than  the  outlined  rank  one  approximation  procedure.  Further,  the  rank  one  procedure 
requires  the  point  source  model  for  the  desired  signal,  which,  in  general,  is  no  longer  the  case  in  radio 
communication  systems. 


5.  OPTIMUM  DETECTION  AND  CROSS-OVER 


It  is  often  the  case  that  the  effect  of  the  spatial  processing  on  subsequent  stages  in  the  receiver 
is  ignored.  In  the  underlying  problem,  it  is  important  to  include  such  an  effect,  specifically,  when  the 
spatial  combiner  varies  from  a  frequency  channel  to  another  for  both  the  suboptimum  and  optimum 
receivers.  The  changes  in  the  power  of  both  the  desired  signal  and  the  jammer,  as  they  marsh  through  the 
spatial  structures  alter  the  sufficient  statistic  for  optimum  detection.  The  sufficient  statistic  for  optimum 
detection  depends  on  the  combiner  response  to  all  sources  impinging  the  aperture,  and  is  given  by 


N 

~  ^  ^m,i 
m=\ 


^*md 


0  ^mj  ^  o  1!! 


(19) 
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where  is  the  combiner  output  of  the  frequency  channel  m  at  symbol  i  after  chip  de-spreading  and 
matched  filtering.  The  scalars  Wn«j(i.e.  w^=w„,".Sd.  Sj  is  the  steering  of  the  desired  in  the  optic  channel) 
and  Wn,j  are  the  responses  of  the  combiner  Wn,  to  the  desired  signal  and  to  the  jammer,  respectively. 


From  (19)  it  is  clear  that  the  response  of  the  combiners  is  required  in  order  to  determine  the 
sufficient  statistics  for  optimum  detection.  Furthermore,  it  may  occur  that  the  spatial  signature  of  the 


jammer  completely  mask  the  desired  signal.  This  is  the  case  when  the  jammer  cross-over  the  DOA  of  the 
desired  source  in  an  optic  channel.  In  consequence,  it  is  necessary  to  keep  control  of  the  adequate  labeling 
at  the  combiners’  outputs.  From  the  three  values  involved  in  each  term  in  (19),  only  the  norm  of  the 
combiner  is  directly  available  to  the  designer;  the  other  two  must  be  estimated  from  data.  It  can  readily  be 
shown  that,  by  taking  the  free  band  in  (15  a)  as  a  reference  band,  the  factor  Fn  for  the  rest  of  the  bands  is 


The  term  corresponding  to  the  spectral  density  of  the  noise  No  can  be  determined  from  the  noise 
eigenvalues  of  .  The  numerator,  in  the  above  equation  can  be  provided  using  the  following  indentity 


P  vv 


(21) 


As  for  the  estimation  of  the  denominator  in  (20),  we  use  two  key  formulations  of  the  design  criterion  for 
the  combiners.  The  first  expression  is  the  MSB  incorporating  (16.b); 


=  (22) 

the  second  relationship,  comes  from  the  definition  of  the  minimization  criterion.  The  MSB  is  lower 
bounded  by  the  residuals  levels  of  jammer  and  noise. 


MSE(n)  >  Jq  Wfjj 


(23) 
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The  equality  holds  for  the  case  where  the  jammer  does  not  shadow  the  desired  signal  and  it  is  properly 
cancelled  in  the  MSE.  Thus,  using  (21),  (22)  and  (23)  in  (20),  we  obtain  an  estimate  of  the  factor  to  be 
used  in  the  decision  variable  for  the  bands  different  from  the  reference  band. 


Note  that  this  factor  is  smaller  than  the  optimum  weigthing  since  the  inequality  (23)  is  used  to 
estimate  the  denominator.  This  factor  accounts  for  the  residual  loss  from  the  spatial  diversity  processing 
due  to  the  presence  of  the  jammer.  In  the  case  where  the  desired  source  is  completely  shadowed  by  the 
jammer,  the  factor  F„  depicts  negligible  contribution  of  the  hit  bands  to  the  decision  vjuiable,  when 
compared  with  the  terms  of  the  bands  free  of  jammer. 


6.  ADAPTIVE  DESIGN 


This  section  describes  an  adaptive  algorithm  to  the  constrained  minimization  described  in  Section 
4.  It  should  be  emphasized  that  the  fundamental  step  is  to  find  the  response  of  the  free-band  combiner  Wf; 
once  this  combiner  is  found  at  each  iteration,  the  update  of  the  rest  of  the  combiners,  including  is 
done  in  a  time  reference  combining  framework,  where  the  output  of  Wf  acts  as  the  reference  waveform. 
As  stated  in  Section  4,  the  presence  of  chip  symbol  modulation  decorrelates  the  jammer  from  one  band  to 
another.  Under  this  premise,  the  band  selection  is  no  longer  crucial  in  the  resulting  performance  of  the 
method. 

The  basic  update  of  the  selected  beamformers  is  done  with  the  instantaneous  gradient  of  the 
Lagrange  objective 
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^  I  (^«  yh,r\  (^n  ^n 

where  (X,,-!)  is  the  Lagrange  multiplier  and 


(25) 


yf,n^y^f,n^f,n 

yh,n^YSi,n^Kn 

are  the  combiners’  output.  The  update  equations  are  then  given  by: 


(26) 


Yif,n^\-'^f,n  Mf,n^f,n^f,n  Kyh,n) 
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The  Lagrange  multiplier  has  to  be  set  such  that  the  new  weights,  prepared  for  snapshot  n+1,  hold  the 
constraint  at  snapshot  n.  In  other  words,  the  updated  combiners  have  to  satisfy  (28). 


The  step-size  parameters  are  set  to  (29), 


(28) 
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The  denominators  in  the  above  equation  are  the  snapshot  powers  defined  by 


(29) 
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WH,n=]Xl^XH,n 


(30) 


Including  the  weight  updates  in  (28)  and  assuming  that  the  past  weights  satisfy  the  constraint,  the 
following  equation  provides  the  selection  for 


+ Afj  (P y^h,n  )(^  (^  2)  —  0 


(31) 
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where 


^h,n  \^h,n—h,n 


(32) 


Since  the  parameter  a  has  to  be  set  equal  to  small  values,  usually  below  0.02,  to  obtain  reasonable  levels 
of  missadjusment,  the  terms  including  a  in  (31)  can  be  neglected  in  favor  of  the  two  terms  which  do  not 
contain  it.  This  results  in 


rn _ 

i^h,n 


(33) 


since  the  constraint  acts  like  an  automatic  control  gain  for  the  system,  the  numerator  of  the  above  formula 
can  be  set  to  a  constant  value.  Note  that  this  constant  will  be,  on  the  average,  the  power  at  the  output  of 
the  free  band.  This  value  was  set  to  one  in  the  simulations  included  in  this  section.  It  is  interesting  to  note 
that,  at  the  steady  state,  as  the  output  powers  tend  to  be  equal,  the  Lagrange  multiplier  will  tend  to  one. 
On  the  other  hand,  at  the  start  up,  when  Pf  and  Ph  are  expected  to  be  different,  the  multiplier  will  be 
smaller  than  one.  In  other  words,  the  Lagrange  multiplier  acts  like  a  flag  indicating  to  the  updating 
algorithm  how  both  output  signals  have  to  be  compared  to  provide  a  proper  error  correction. 

A  realistic  set-up  of  the  algorithm  demands  some  smoothing  of  the  instantaneous  powers  defined 
previously.  The  resulting  performance,  in  terms  of  convergence  rate,  does  not  show  high  sensitivity  to 
this  parameter,  providing  it  is  above  a  threshold.  For  different  scenarios,  the  values  of  the  smoothing 
factor  which  are  above  6.8  shown  small  impact  in  the  learning  curve. 

In  Figure  2  depicts  the  learning  curve  (i.e.  the  smoothed  MSB  versus  number  of  snapshots)  for 
two  different  scenarios,  which  differ  in  the  angular  separation  between  the  desired  source  and  the  jammer. 
The  rest  of  parameters  are:  FDSS  of  6  bands,  T)  equal  to  50%,  5  sensors  half  wavelength  ULA  array,  the 
desired  and  jammer  signals  are  1  dB  and  10  dB  respectively,  a=0.01  and  P=0.9.  From  this  Figure,  it  is 
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clear  that  the  adaptive  algorithm  exhibits  good  behavior  which  makes  it  a  candidate  alternative  to  frame 
processing  method  suitable  for  TDMA  which  has  been  reported  in  Section  4. 


7.  COMMENTS  AND  EXTENSSIONS  OF  THE  PROCEDURE 


It  is  important  to  note  that  the  proposed  method  is  based  merely  on  the  coded  stracture  of  the 
desired  signal  in  the  frequency  domain.  We  should  stress,  however,  that  the  applicability  of  this  method  is 
granted  irregardless  of  the  nature  of  the  code,  i.e.,  wether  it  is  being  pure  repetitive  or  otherwise.  With  the 
same  perspective,  it  is  easy  to  realice  that,  as  long  as  the  desired  signal  has  a  stracture  (  we  refer  to  “a 
stracture”  as  any  repetition  in  any  diversity  axis  like  encode,  frequency,  time,  etc.),  optimum  spatial 
processing  is  attainable. 

From  the  above  perspective,  spatial  processing  methods,  referred  to  as  blind,  are  not  enterely 
blind,  since  some  information  concerning  the  signal  stracture  is  used  at  the  receiver  to  obtain  optimum 
processing  .  In  communications  systems,  this  stracture  is  provided  for  proper  matching  between 
modulation  and  demodulation  schemes  at  both  ends  of  the  communications  channels. 

A  way  around  to  describe  the  underlaying  advantages  of  FD  systems  for  blind  beamforming  is  to 
traduce  the  signal  structure  to  the  property  of  exact  prediction.  The  stracture  of  the  FDSS  transmitted 
signal,  in  essence,  has  allowed  the  desired  signal  to  be  exactly  predicted  from  one  band  to  another,  and 
thereby  serves  as  a  reference.  Since  this  scheme  is  using  predictability  in  the  frequency  domain  it  is 
evident  that  some  correspondence  has  to  exist  in  the  time  domain.  This  possibility  was  outlined  in  Section 
2,  when  it  was  proposed  to  divide  the  chip  signal  in  several  segments  within  a  symbol  and  to  transmit 
every  slot  in  a  different  frequency  band.  The  time  domain  dual  is  the  transmissions  of  these  slots 
sequentially  in  time;  time  diversity,  to  counteract  partial  time  jamming  [7].  With  time  diversity,  the 
segments,  within  an  information  symbol,  can  be  predicted  which  are  then  used  to  form  the  blind 
combiner.  For  clarity  of  the  presentation,  let  us  assume  that  the  chip  signal  within  a  symbol  is  divided  in 
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two  segments,  as  indicated  in  Figure  3.  Note  that  the  pure  repetition  code  is  also  included  within  this 
general  scheme  (i.e.  the  same  chip  code  is  used  twice  per  symbol  interval). 

Since  the  described  structure  in  the  time  domain  is  exactly  the  same  as  that  of  the  FDSS  in  the 
frequency  domain,  the  optimum  spatial  combiner  for  the  corresponding  desired  signal  can  be  derived 
from  the  minimization  of  ^  defined  as  (34), 


=  £: 
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where  2Li.i  tind  2£n.2  are  the  snapshots  corresponding  to  segments  one  and  two  respectively.  This 
formulation  is  adequate  for  partial  time  jamming,  as  the  desired  signal  temporal  structure  allows  exact 
prediction.  More  interesting  is  the  case  where  the  minimization  of  (34)  is  carried  over  a  single  combiner, 
which  is  possible  if  the  jammer  remains  uncorrelated  from  time  segment  to  segment , 


^n,2 


(35) 


Since  the  jammer  and  the  noise  are  assumed  uncorrelated  for  the  two  segments  selected,  the 
optimum  combiner  will  select  only  the  signal  which  is  exactly  predicted  by  a(n),  i.e.  the  desired  signal.  Of 
course,  the  additional  constraint  used  for  FDSS  has  to  be  added  in  order  to  avoid  the  trivial  solution  to  the 
minimization  of  (35).  The  exact  prediction  required  for  the  procedure,  outlined  previously,  can  be  also 
provided  by  the  channel  encoder.  Note  that  using  the  same  combiner  precludes  optimal  detection,  but 
reduces  the  complexity  in  the  suboptimum  detector. 


For  further  illustration  of  the  notion  of  signal  prediction  and  diversity  we  selected  below  the  case 
of  an  unmodulated  carrier  at  baseband  frequency  f©.  Assuming  that  the  maximum  time  for  coherence  in 
the  jammers  is  smaller  than  to  (  This  is  the  case  when  to  is  greater  than  the  inverse  of  the  interferer’s 
symbol  rate),  the  minimization  objective  is  (36). 
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(36) 
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This  objective,  which  is  constrained  by  the  non-zero  value  of  the  cross-correlation  between  the 
two  outputs  produces  adequate  nulling  of  the  jammer  at  the  combiners’  output.  Defining 


R  =  )£  =  Ek-nK.”r,+to 

the  optimum  combiner  is  the  eigenvector  associated  with  the  minimum  eigenvalue  of 


(37) 


Ryv  =  APw  (38) 

Figure  4  shows  an  example  using  the  exact  prediction  property  of  an  unmodulated  carrier.  The 
top  part  of  this  Figure  shows  the  beamformer  response  (top)  whereas  the  bottom  part  depicts  the 
beamforming  scheme,  which  reflects  the  low  complexity  of  the  reference  loop.  The  scenario  contains  two 
BPSK  interferers  (10  and  3  dB  of  SNR)  at  symbol  rate  greater  than  to’’  with  DOAs  equal  to  0  and  10°, 
together  with  the  desired  unmodulated  carrier  (0  dB.)  impinging  from  -20°.  The  baseband  frequency  is 
0.25  for  all  sources.  The  adequate  response  of  the  combiner  designed  from  500  snapshots  (250  each  set), 
using  the  afforementioned  procedure  is  clearly  evident. 

It  is  noteworthy  that  this  system  requires  only  the  nominal  frequency  of  the  desired  signal  and 
does  not  need  phase  synchronism  which  is  very  convenient  in  practice.  This  allows  a  great  saving  in 
complexity  and  threshold  margin  for  very  low  SNRIs. 

Finally,  note  that  since  frequency  structure  effectively  remove  partial  band  jamming,  time 
structure  does  the  same  with  partial  time  jamming.  In  consequence,  if  the  structure  of  the  transmitted 
signal  has  simultaneously  time  and  frequency  diversity,  the  receiver  will  be  robust  against  these  two 
undesired  jamming  effects. 
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8.  SIMULATIONS 


The  most  important  issue  concerning  the  usefulness  of  spatial  combining  for  frequency  diversity 
is  the  adequate  weighting  of  every  term  of  the  sum  forming  the  decision  variable  Am  in  (19).  This  factor 
reflects  the  limitations  of  the  combiner  when  the  jammer  shadows  globally  or  partially  some  bands  of  the 
FD  system.  The  estimate  of  this  factor  is  proposed  in  (23),  as  is  shown  in  Figure  5  presents  this  estimate 
for  different  DOA  separation  between  the  jammer  and  the  desired  signal.  For  every  separation,  a  set  of 
400  snapshots  has  been  used.  The  variation  of  the  DOA  separation  was  set  from  0"  (complete  shadowing) 
up  to  20‘’.The  SNR  of  the  jammer  was  in  all  cases  10  dB  (SJR=-10dB).  No  chip  modulation  has  been  used 
in  the  simulations;  in  consequence,  the  jammer  is  fully  correlated  from  band  to  band,  and  the  performance 
reported  can  be  considered,  in  practice,  as  a  lower  bound.  The  desired  signal  is  BPSK  modulated  and  its 
power  is  changed  from  10  dB  down  to  0  dB  (step  1  dB).  The  desired  signal  permanent  location  was  set 
equal  to  10°  from  the  broadside.  The  FDSS  system  has  six  frequency  bands  out  of  which  bands  1,  2  and  4 
are  hit  by  the  jammer  (T|=  50%)  at  10  dB  of  Jammer  to  Noise  Ratio  (JNR)  each.  The  labeling  of  the  free 
and  hit  bands  was  done  from  a  power  detector.  The  aperture  was  an  ULA  array  of  5  sensors.  The  estimate 
of  the  shadowing  factor  Fn  improves  when  the  SNR  of  the  desired  is  high,  that  is  Fn  is  close  to  zero  for 
zero  degrees  of  sources  separation  and  close  to  one  when  the  jammer  is  impinging  on  the  aperture  from 
secondary  lobes  (for  this  case  the  spatial  bandwidth  of  the  aperture  at  the  central  frequency  was  20  °). 
Performance  deterioration  for  low  signal  to  noise  ratio  is  expected,  since  the  difference  between  the  free 
and  the  hit  bands  decreases,  with  high  levels  of  noise.  When  factor  F„  was  used  in  the  detector  statistics, 
the  Bit  Error  Rate  (BER)  was  always  below  or  equal  (complete  shadowing)  to  the  resulting  one  for  the 
suboptimum  detector  (BER  corresponding  to  SNRd  plus  the  number  of  free  bands  (lO.logS  dB),  plus  the 
number  of  spatial  diversity  channels,  (lOlogS  dB)). 

It  is  also  interesting  to  measure  for  the  same  scenario  the  change  in  the  SNR  over  the  different 
bands  at  the  combiners’  outputs  for  different  SJRs.  Figure  6  shows  (top)  for  the  hit  bands  and  (bottom) 
for  the  free  bands  the  average  of  the  actual  SNR  for  each  group,  together  with  the  impact  of  the 
corresponding  shadowing.  The  input  SNRd  is  set  to  1  dB.  It  is  important  to  note  that  the  contribution  of 
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the  SNR  in  the  free  bands  increases  with  increased  jammer  power  and  with  the  shadowing  effect.  The 
worst  case  is  represented  by  a  30  dB  jammer  above  the  desired  in  every  hit  band  and  DOA  equal  to  the 
same  as  desired  signal.  The  result  is  1  dB  less  than  the  case  when  there  is  two  degrees  of  separation 
between  the  jammer  and  the  desired  signal.  This  plot  provides  evidence  that  the  new  front-end  is  able  to 
resort  at  any  time  to  the  suboptimum  detector.  This  is  because  the  free  bands  remain  almost  unaffected  by 
the  spatial  combiner  processing.  The  role  of  spatial  processing  for  improving  SNR  in  the  hit  bands  should 
be  well  appreciated.  This  traduces  in  an  effective  reduction  in  factor  T|  which,  for  this  aperture,  reduces  to 
zero  when  the  separation  of  the  desired  and  the  source  is  greater  than  6  degrees. 

9.  CONCLUSIONS 


It  has  been  shown  that  the  use  of  diversity,  basically  frequency  diversity,  in  the  sense  of  several 
replicas  of  the  desired  signal,  allows  the  optimum  design  of  the  spatial  combiners  blind  to  the  source 
waveform.  In  a  wide  sense,  it  can  be  stated  that  whenever  the  desired  signal  along  any  diversity  axis  has 
an  exact  representation  over  diversity  components,  the  spatial  combiner  can  be  optimally  designed 
without  additional  diversity  slots  for  a  reference  signal.  The  propossed  design  procedure  minimizes  the 
MSE  between  two  selected  diversity  components,  as  in  an  exact  prediction  problem,  with  the  constraint  of 
a  cross-correlation  factor  different  from  zero.  It  is  remarked  that  the  diversity  could  be  time,  frequency, 
encode  or  any  mixture  of  them.  The  procedure  is  blind  to  the  desired  signal  waveform.  Since  the  spatial 
diversity  may  collapse,  depending  if  the  jammer  may  sign  in  all  diversity  components,  the  optimum 
detector  demands  for  an  estimate  of  this  effect  in  order  to  properly  adjust  the  likelihood  in  every  diversity 
branch.  To  cope  with  this  problem,  an  estimate  has  been  derived  from  the  design  parameters.  This 
estimate  almost  loads  the  detector  to  operate  only  on  the  jamming  free  diversity  branches  when  the 
jammer  shows  strong  presence.  This  is;,  in  essence,  an  automatic  convergence  to  the  so-called  sub¬ 
optimum  receiver.  It  has  been  outlined  that,  whenever  any  kind  diversity  allows  exact  prediction  of  only 
the  desired  signal  over  two  diversity  components,  blind  spatial  processing  is  always  available. 
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FIGURE  CAPTIONS 


Figure  1.  General  transmitter  scheme  for  spread  spectrum  digital  communications. 

Figure  2.  Learning  curve  for  the  adaptive  algorithm,  i.e.  the  smoothed  MSB  versus  snapshots  used,  for 
DOA  separations  of  1°  (top)  and  15®  (bottom)  between  the  desired  signal  and  the  jammer.  The  scenario 
and  algorithm  parameters  are:  6  bands,  50%  bands  jammed,  IdB  and  lOdB  of  SNR  for  desired  signal  and 
jammer  respectively,  a=  0.01,  P=0.9  and  5  sensors  half  wavelength  ULA  array. 

Figure  3.  The  code  stracture  within  time  diversity:  A  segment  of  the  desired  signal  can  be  exactly 
predicted  from  the  previous  segment. 

Figure  4.  Beamforming  response  based  on  the  exact  prediction  of  the  desired  signal.  Desired 
(unmodulated  carrier)  and  jammers  are  labeled  with  lines  in  both  plots  on  top.  Bottom  the  beamforming 
architecture  with  the  reference  loop. 

Figure  5.  The  shadowing  factor  Fn  as  a  function  of  the  DOA  separation  of  the  desired  versus  the  jammer 
for  several  SNR  (SNRd)  of  the  desired  signal. 

Figure  7.  Average  actual  SNR  for  the  hit  bands  (top)  and  for  the  free  bands  (bottom)  versus  Signal  to 
Jammer  Ratio  (SJR),  for  different  DOA  separation  between  the  desired  source  and  the  jammer. 
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Abstract 


This  paper  proposes  the  subband  adaptive  array  processing  method,  which  is  used  to  mitigate  both 
inter-symbol  interference  (ISI)  and  cochannel  interference  (CCI)  in  digital  mobile  communications. 
The  analysis  filters  enhance  the  correation  of  multipath  rays  in  each  subband.  This  enhancement  is 
blind  in  the  sense  that  no  a  priori  knowledge  of  the  arriving  signals  is  required.  With  the  increased 
coherence,  the  desired  signal  can  be  effectively  equalized  by  subsequent  spatial  processing.  Further, 
the  CCI  signals  and  their  multipaths  can  be  suppressed  with  fewer  array  degrees-of-freedom.  The 
effect  of  QMF  and  DFT  filter  banks  on  multipath  correlation  is  delineated. 
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1.  Introduction 


In  land  mobile  communications,  a  signal  transmitted  by  a  mobile  station  arrives  at  a  base  station 
reflected  and  scattered  by  surroundings,  which  causes  multipath  fading.  As  the  demand  is  increasing 
for  multimedia  communications,  mobile  communication  is  developing  towards  high-speed  digital 
wireless  networks,  where  the  communication  channels  is  frequency-selective,  and  the  inter-symbol 
interference  (ISI)  becomes  highly  pronounced.  Another  important  problem  in  mobile  communication 
is  cochannel  interference  (CCI),  which  is  generated  due  to  the  frequency  reuse  in  cellular  systems. 

Adaptive  array  implementing  spatial  or  spatial-temporal  equalizations  or  diversity  techniques 
prove  useful  in  suppressing  both  of  the  ISl  and  CCI,  leading  to  increased  capacity  and  range  [1-8]. 
The  spatial-temporal  equalization  can  be  achieved  by  space-time  adaptive  processing  (STAP)  to 
effectively  mitigate  ISI  and  CCI.  Such  a  scheme  is  composed  of  an  integrated  adaptive  array  and 
temporal  equalizer  to  perform  optimum  jointly  spatial  and  temporal  signal  processing.  However, 
solving  both  the  CCI  and  ISI  problems  simultaneously  by  conventional  STAP  methods  is  yet 
difficult.  These  methods  require  either  large  scale  matrix  inversion  or  recursive  calculation,  or  a 
cascaded  form  of  CCI  and  ISI  cancellers  [6-8]. 

In  this  paper,  we  propose  an  efficient  subband  adaptive  array  processing  method  that  utilizes  filter 
banks  to  mitigate  both  the  CCI  and  ISI  effects  in  land  mobile  communications.  A  subband  adaptive 
array  has,  in  effect,  the  same  function  as  an  STAP  system,  while  the  implementation  is  much  easier. 
The  subband  adaptive  array  provides  sub-optimal  performance  in  the  context  of  frequency-spatial 
signal  processing,  which  enhances  the  signal  correlation  between  multipath  rays  prior  to  processing. 
In  subband  adaptive  arrays,  the  frequency  band  of  the  received  signals  is  divided  into  smaller  bands 
by  the  use  of  filter  banks  [9].  Proper  analysis  filters  allow  a  significant  increase  in  the  signal 
correlation  between  the  multipath  rays  within  each  subband.  Note  that  such  an  increase  is  blind  in  the 
sense  that  it  does  not  require  a  priori  knowledge  of  the  arriving  signals.  As  a  result  of  increased 
correlation,  the  multipath  fading  associated  with  both  the  desired  and  the  interference  signals  is 
reduced,  thus  yielding  faster  convergence  of  the  adaptive  weight  vector  as  well  as  fewer  degrees-of- 
freedom  (DOF’s)  required  for  adequate  equalization  [10,1 1]. 

So  far,  subband  adaptive  array  methods  have  been  proposed  for  wideband  signal  processing. 
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Subband  adaptive  arrays  have  the  advantages  of  reducing  the  processing  rate  [14]  as  well  as 
achieving  rapid  convergence  [12,  15].  Other  kind  of  coherent  subspace  transformation  methods  [13, 

1 6]  also  demonstrate  these  two  advantages. 

In  this  paper,  the  discussion  is  mainly  focused  on  the  signal  correlation  enhancement  in  a 
frequency-selective  fading  environment,  rather  than  the  conventional  consideration  on  subband 
signal  processing  to  reduce  the  signal  bandwidth  and  computational  complexity,  although  such 
advantages  are  also  shared  by  our  proposed  method.  The  signal  correlation  enhancement,  in  essence, 
converts  a  frequency-selective  fading  problem  into  a  flat  fading  one.  In  this  case,  the  ISl  problem  is 
highly  mitigated  as  the  faded  signal  can  be  easily  equalized  by  a  subsequent  spatial  combiner.  Since 
the  CCI  components  in  each  subband  are  easily  removed  by  the  respective  narrowband  adaptive 
array,  subband  adaptive  array  techniques  can  simply  solve  both  CCI  and  ISI  problems. 

Although  many  transformation  methods  can  convert  a  wideband  signal  into  multi-channel 
narrowband  signals,  it  is  shown  in  this  paper  that  different  methods  provide  different  signal 
correlation  enhancement  performance.  The  importance  to  properly  choose  the  analysis  and  synthesis 
filters  to  obtain  good  equalization  performance  by  subband  signal  processing  is  emphasized. 

The  paper  is  organized  as  follows.  In  Section  II,  we  first  present  the  framework  of  -adaptive  array 
processing  in  the  presence  of  multipath  rays.  It  is  shown  that  multipath  rays  of  a  signal  can  be  treated 
as  a  single  ray  when  their  cross-correlation  is  relatively  high.  The  signal  correlation  enhancements  of 
multipath  signals  by  using  different  filter  banks  are  compared  in  Section  III.  In  Section  IV,  the 
analysis  of  the  subband  array  performance  is  given.  In  Section  V,  simulation  examples  are  provided 
to  confirm  the  effectiveness  of  the  proposed  method. 

II.  Adaptive  Array  Framework  in  the  Presence  of  Multipath  Rays 

Without  a  loss,  of  generality,  we  consider  the  simple  case  where  a  desired  signal  and  a  CCI,  each 
with  two  rays,  are  incident  on  an  array  of  AT  sensors.  The  received  array  vector  is  described  as 

x(/)  =  A(©;,  )s  o  (0  +  A(©  /  )s  /  (/)  +  n(/) ,  ( 1 ) 

where  A(©c)  =  [a(0^,),a(0y,,)]  is  the  collection  of  the  steering  vectors  for  the  desired  signal,  and 
S/)(0  =  [^di(0>'Sd2(0]^  is  the  vector  which  represents  the  associated  complex  envelopes  at  the 
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presumed  phase  reference  point  of  the  array.  Similarly,  A(©/)  =  [a(6’;,),a(^/j)]  and 

S/(0  =  U/i(0.'^y2(0]^are  the  corresponding  vectors  for  the  CCI  signal.  It  is  assumed  that  the  desired 
signal  and  the  CCI  signal  are  independent.  We  also  assume  that  the  elements  of  the  noise  vector 
n(/)  are  i.i.d.  Gaussian  random  processes  with  zero-mean  and  covariance  matrix  cr^I ,  where  I  is 
the  A^x  //identity  matrix.  is  used  as  the  direct  ray  and  considered  the  signal  of  interest  (SOI). 

The  multipath  ray  5^,  (/)  can  be  separated  into  two  components,  one  is  coherent  with  the  SOI  and 
the  other  is  orthogonal  to  the  SOI  [17],  i.e., 

~  ^[>(.0  (2)  * 

where  ^  (O^^i  (0]  =  0 .  2nd 

Pd  ~  wy 

D1 

is  the  signal  correlation  between  the  two  rays  whose  power  are  £rc(=o’|ii)  2nd  ,  respectively, 

-^Dz^^Dx-  denotes  statistical  average  and  (•)*  denotes  complex  conjugate.  The  signal 
covariance  matrix  is  defined  as 


R'„  =  £Is„(/)s“(()] 


Pd^d 


p'd^d  \Pi^dV 


] 


(4) 


where  denotes  complex  conjugate  transpose.  Accordingly,  the  spatial  covariance  matrix  is  given  - 
by 


Rr>y>  =  A(0y,)£[s^(Os2(O]A''(©^) 


1 

Pd^d 


=  o-fl[a(i9r„)  +  PcAD^i^Dz)^  +  0  “IPof 


(5) 

The  first  term  of  the  right  side  of  (5)  is  constructed  from  combined  signal  vector 
Sc (0  =  [a(^ci )  p'd^d^^^dz (0  .  where  the  quantity  in  the  brackets 
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(6) 


represents  the  generalized  steering  vector  of  SOI.  The  second  term  in  (5)  is  the  contribution  to  the  ISI 
from  the  orthogonal  components  s^{t) .  This  contribution  vanishes  asymptotically  as  \p^\ 
approaches  1 . 

By  assuming  that  the  reference  signal  r{t)  is  an  ideal  replica  of  the  SOI  then  the 

correlation  vector  between  the  data  vector  x(/)  and  the  reference  signal  rif)  is 


=  £[x(/)/'*  (0]  =  (01  = 

which  coincides  with  the  generalized  steering  vector  of  the  SOI,  given  in  (6).  Thus,  the  minimum 
mean  square  error  (MMSE)  optimum  weight  vector  is  given  by  [18] 

=  =gR''ao  (8) 

where  g  =  o-o(l-o-^a^R^a^)  is  a  constant, 

R^,  =  ^[xCOx''  (/)]  =  R;,  +0-’ a^a^  (9) 

RAf=(l~|Pcl  (^02)  ^/,V  ■' 

R;v  is  the  equivalent  covariance  matrix  of  all  interfering  components,  and 

R,v  =R;,  +£[n(0n''(0]  =  R;,  ■  (H) 

R„  is  the  covariance  matrix  of  the  CCI  signal  vector,  which  is  defined  in  (14).  By  using  the  matrix 
inversion  formula,  we  obtain 

Ra'  ~  Rw  ~  ^R/iV®(^£I2)®  (^D2)R/,V 

where 

_ 0~|Pz)|  )^Di _  ^12^ 

l+«(l~|Po|  )<7d2®  (^i)2)R/.V^(^r)2) 

which  monotonously  decreases  as  |p^|  increases. 


From  the  above  analysis,  the  array  weight  vector  will  be  less  influenced  by  the  orthogonal 
component,  5^(0 .  as  |/j^|  increases.  For  the  case  of  \p^\  reaching  a  unit  value,  the  two  multipath 
rays  can  be  considered  a  single  ray,  and  they  can  be  handled  by  a  single  DOF.  Consequently,  the 
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frequency-selective  fading  problem  is  converted  into  a  flat  fading  problem.  Thus,  the  channel  is 
equalized  by  array  processing  and  the  ISI  problem  is  solved. 

All  CCI  components  are  to  be  suppressed  by  the  spatial  signal  processing.  Therefore,  these 
components  can  be  represented  by  using  eigen-decomposition  methods.  The  covariance  matrix  of  the 
CCI  signal  vector,  denoted  by  R/, ,  is  expressed  as. 


R,,  =A(e,)£[s,(/)sf  =  (14) 

/si 

where  Aj.(i=\,2)  are  the  eigenvalues  of  and  u^,.  are  the  associated  eigenvectors.  The  two 
eigenvalues  are  given  by  [19] 


^11.2  =  yI®"/!  +cr/2  Re(/7;/?;)] 


\cr i\  +o’/2  ■*‘2o'/iC7/2  Re(/3y^y)] 


] 


(15) 


where  a]^  and  cr^,  are  the  powers  of  the  two  rays  of  the  CCI  signal,  respectively,  p,  is  the  signal 
temporal  correlation  coefficient  between  the  two  rays,  and  p,  is  the  spatial  correlation  coefficient  of 
the  rays,  defined  as 


a^(g„)a(g,2) 

||a(^/.)||||a(<9/3)|| 


(16) 


It  is  clear  that  becomes  0  when  either  or  both  of  \p,\  and  \p,\  become  1.  For  these  cases,  the 
two  rays  behave  as  a  single  ray.  Therefore,  while  the  array  often  needs  two  DOF's  to  effectively 
suppress  the  two  rays  with  low  signal  correlation,  it  suffices  only  to  use  a  single  DOF  when  the  two 
rays  are  highly  correlated  statistically  and/or  spatially. 


III.  Signal  Correlation  Enhancement  by  Using  Filter  Banks 

In  the  previous  section,  we  discussed  the  adaptive  array  framework  in  the  presence  of  multipath 
rays.  It  is  shown  that,  the  increase  of  the  correlation  between  signal  multipath  rays  converts  the 
frequency-selective  problem  into  a  flat  fading  problem,  so  that  the  signal  is  equalized  by  array 
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processing.  When  this  correlation  is  very  high,  the  multipath  rays  can  be  treated  as  a  single 
equivalent  ray.  A  generalized  beam  is  formed  for  the  multipath  rays  of  the  desired  signal  by  a  single 
array  DOF  so  that  the  rays  are  combined  in  the  manner  of  maximum  ratio,  and  the  multipath  rays  of 
each  CCI  signal  will  be  suppressed  by  a  single  array  DOF. 

This  section  discusses  the  signal  correlation  enhancement  of  multipath  rays  by  subband  signal 
processing.  We  first  briefly  review  the  concept  of  subband  signal  processing  by  using  filter  banks, 
then  give  the  basic  principle  of  coherency  enhancement  via  subband  partition  using  ideal  filters. 
Next,  the  implementation  based  on  DFT  and  QMF  filter  banks  is  discussed  and  the  performances  of 
signal  correlation  enhancement  are  compared. 

A.  Subband  signal  processing  by  using  filter  banks 


A  filter  bank  is  a  set  of  filters  linked  by  sampling  operators  and  sometimes  by  delays.  Fig.l  shows 
the  block  diagram  of  a  filter  bank.  The  r-transform  of  the  output  y{n)  is  expressed  as 


y(z)  =  nz)Xiz)  +  Y,A,{z)X{e^-’’'‘'^z) 

(17) 

where  X{z)  is  the  z-transform  of  ;c(w), 

1  M 

M  *=i 

(18) 

M  jtsi 

(19) 

and  M  is  the  number  of  subbands,  and  we  assume  the  decimation  rate  is  the  same.  This  filter  bank  is 

called  maximally  decimated  filter  bank.  The  first  term  of  the  right  side  of  (17)  is  the  non-aliasing 
component,  and  the  last  term  is  the  aliasing  component.  To  realize  the  perfect  reconstruction  (PR)  of 

the  filter  bank,  it  is  required  that 
# 

T{z)  =  r"'" 

(20) 

and 

At(z)  =  0  {ork=  1,2, 

(21) 

where  an  arbitrary  integer. 
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Under  the  PR  condition,  the  output  of  the  filter  bank  is  a  delayed  version  of  the  input  of  the  filter 
bank  without  distortion,  or 

yin)  =  x(n-no).  (22) 

If  a  filter  bank  consists  of  ideal  filters  is  used,  that  is 

’■  ’  elsewhere 

then  clearly  the  filter  bank  satisfies  the  PR  condition.  Although  such  a  filter  bank  is  not  casual,  it 
provides  a  very  simple  expression.  The  expression  will  be  used  to  analyze  the  signal  correlation 
enhancement  in  sub-section  B. 

A  useful  type  of  PR  filter  bank  is  DPT  filter  bank,  which  is  based  on  DFT  and  IDFT  transforms,  as 
is  shown  in  Fig.  2.  The  DFT  filter  bank  is  discussed  in  sub-section  C. 

Quadrature  mirror  filter  (QMF)  banks,  or  more  practically,  pseudo-QMF  banks,  are  often  used  to 
realize  efficient  subband  signal  processing  with  perfect  or  near-perfect  reconstruction  (NPR)  [9,  20, 
21].  A  QMF  filter  has  the  same  construction  as  Fig.  1,  where  in  a  QMF  filter  bank  each  of  the 
analysis  or  synthesis  filters  can  be  realized  by  an  FIR  filter  with  finite  taps.  In  a  pseudo-QMF  bank, 
the  filter  bank  channels  are  formed  by  a  series  of  equidistant  frequency  shifts  of  an  appropriate 
prototype  filter,  and  the  transfer  functions  of  adjacent  channels  are  approximately  power 
complementary  between  their  center  frequencies.  We  will  discuss  the  pseudo-QMF  banks  in  sub¬ 
section  D. 


B.  Basic  principle  of  coherency  enhancement  via  suhhand  partition 


Consider  a  band-limited  stochastic  signal  s(t) .  For  the  simplicity  of  analysis,  we  assume  the 
power  spectrum  density  function  (PSDF),  denoted  by  p{f) ,  is  flat  within  bandwidth  B , 


Pif) 


elsewhere. 


(24) 


The  associated  autocorrelation  function,  denoted  by  r(r) ,  is  obtained  by  the  Fourier  transformation 
of  />(/),  as 
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TIT 

If  we  equally  partition  p{f)  into  A/ subbands,  such  as 

p(/)  =  Ep'*^(/), 

k^O 

where 


{.  B  k  o 

1, - +  B<  f  < - + 

0  2  A/  2 


B  k  +  \ 


-B 
M  » 


elsewhere 

then  the  autocorrelation  function  for  Ath  subband  signal  will  be 

sin(^^)  a 

r**>  (r)  =  £{5“>  (/  -  r)}  = - . 

TtT 

The  signal  correlation  coefficient  of  s{t)  and  its  delayed  replica  s{t  -  r)  is  then 
^  ^  r(0)  nBr 


(25) 


(26) 


(27) 


(28) 


(29) 


where  a  (|al=l)  is  a  factor  used  to  express  the  propagation  phase  difference  between  ^(0  and 
s{t  -  r) .  Similarly,  for  the  kth  subband  signal  we  have 


p<*)(^)  =  .  (30) 

^  r<*\0)  tcBtIM  ^  ’ 

By  comparing  (29)  and  (30),  it  becomes  apparent  that  the  signal  correlation  of  the  subband  signal  is 
widened  A/ times  in  the  time  scale.  We  note  that  |p^*’(r)|  does  not  change  with  A,  which  implies  that 

the  signal  correlation  is  equally  enhanced  for  the  different  subbands.  For  example,  when  r  =  1/5  we 
have  /j(l/5)  =  0,  which  means  two  signals  are  uncorrelated.  When  the  signal  is  divided  into  M 

subbands,  then  (1/5)1  =  which  equals  to  0.9936  for  A/ =16.  Therefore,  the  two 

'  ‘  nIM 

uncorrelated  rays  become  approximately  coherent  in  the  subbands.  The  correlation  |/j'*^(r)|  is 
plotted  versus  Afand  is  shown  by  the  solid  line  in  Fig.  3  for  5r  =  1 . 
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C.  Subband  partition  via  DFT  filter  banks 


In  the  above  sub-section,  we  considered  the  filter  bank  based  on  ideal  filters.  However,  for  causal 
systems,  the  transfer  function  of  an  FIR  filter  does  not  have  the  shape  of  a  rectangular  pulse.  The 
most  practically  used  filter  banks  are  based  on  DFT  and  QMF  filter  banks  to  realize  perfect 
reconstruction  or  near-perfect  reconstruction  [9,  21],  In  this  sub-section  we  analyze  the  performance 
of  signal  correlation  enhancement  when  DFT  filter  banks  are  used. 

The  discrete  Fourier  transformation  (DFT)  of  a  data  sequence  x{n)  is  defined  as 

Xik)  =  Yx(r7W'"  ,  k^0X-,N-l  (31) 

«=0 

where  W  =  .  The  Ath  subband  output  of  the  DFT  is 

_y<*)(/)  =  *;^'jc(r  -  ^  0  <  A  <  M  - 1  (32) 

«=0 

where  is  the  sampling  duration.  Here  we  assume  critical  sampling,  so  that  =  l/B .  The  impulse 
responses  of  Ath  subband  DFT  analysis  and  synthesis  filter  can  be  described  as 


0<n<M-l,  0<k^M-\. 

Then,  the  z-domain  transfer  functions  of  the  filters  are  described 

.  »=0 

0<A<M-1 

respectively,  and  the  frequency  response  of  the  Ath  analysis  filter  is 

//<*)(/)  =  ^'g-;2»'/g-y2;r(.W-l-;0*/W 
«=0 

_  sin(;r M{f  —  kl M)) 

where  /  denotes  the  normalized  frequency  with  respect  to  the  signal  bandwidth  B. 
From  (35)  and 


(33) 


(34) 


(35) 
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(36) 


the  autocorrelation  function  and  signal  correlation  coefficient  of  subband  filter  output  signal  is 
obtained  as 

-* 

.^nsO  msO 

Af-1A/-I  fj  _  fyt 

=  Z 

lieO  maO  ^ 

and 

/?<*>(r)  =  ar‘*>(r)/r<*>(0).  (38) 

The  signal  correlation  coefficient  is  plotted  in  Fig.4  for  the  different  subbands.  It  is  seen  that  the 
curve  is  flat  only  when  r/r^  is  an  integer.  At  the  other  cases,  the  signal  correlation  coefficient  has 

M 

the  largest  value  at  jfc=0,  and  the  coefficient  is  very  low  at  i=  — ,  particularly  when  the  modulo  of 
rjTs  is  0.5. 

The  correlation  coefficient  bet%veen  the  two  rays  is  given  as  the  dotted  line  in  Fig.'  3  for  Bt  =  \ 
and  k  =0.  |p^*^(r)|  is  approximately  0.935  when  il/=16.  From  Fig.  3  it  is  clear  that  the  DFT  filter 

banks  provide  meager  signal  correlation  enhancement.  From  Fig.4  it  is  seen  that  |/7^*'(r)|  is  even 

lower  when  r/T^  is  not  an  integer  and  k  is  close  to  M2.  Thus,  filter  banks  with  more  powerful 
correlation  enhancement  are  preferable  to  reduce  the  number  of  subband  partitions. 

D.  Subband  partition  via  QMF  filter  banks 

We  now  discuss  the  signal  correlation  enhancement  using  pseudo-QMF  filter  banks.  The  impulse 

♦ 

responses  of  the  pseudo-QMF  analysis  and  synthesis  filters  hf^-iii)  and  /^.(w)  respectively,  are 
cosine-modulated  versions  of  a  real  and  linear-phase  impulse  response  of  the  prototype  filter  h{ri) , 
i.e. 


67 


(39) 


h,.  (n)  =  2Kn)  cos((2;r‘+l)  ^  ^ 

/,.(«)  =  2Mn)  cos((2*’+l)  (n  -  -  (-!)*• 

^  J.  4 

0<«  <//-!,  0<k'<M'-\, 


K-\ 

where  N  is  the  number  of  FIR  taps.  Denote  H^z)  -  ^/:(«)z‘''  as  the  z-transform  of  h{n) ,  the  z- 

/;=0 

transforms  of  h^.{n)  are 


(z)  =  a,.c^.H{zW  ^  )  +  a].c].H{zW  ^  ) 


(40) 


where  W  =  e'^'' ,  a^.-  e’^‘' ,  c*.  =  W 


t. 

■  ^  ,  and  $1,.  =  (-1)  — .  Due  to  the  mirror  effect,  the 
4 


signal  correlation  is  very  low  for  a  large  k’  and  non-integer  tJT^.  Therefore,  direct  utilization  of 
QMF  subband  filter  bank  does  not  bring  adequate  signal  correlation  enhancement  for  different 
subbands.  To  solve  this  problem,  we  partition  each  subband  analysis  filter  into  two  sub-subband 
filters,  //^*  ''(z)  and  H^^'^\z).  Consequently,  in  our  modified  QMF  subband  processing  scheme, 
the  analysis  filter  bank  consists  of  lAf  subband  coding  (SBC)  filter  banks  while  the  synthesis  filter 
bank  consists  of  Af  QMF  banks. 

To  simplify  and  unify  the  notation,  we  let  M=2M\  and  denote  the  subband  index  {k',  i)  into  k  so 
that 


_  (  it'  /  = 

“  lit'+A/’  /  = 


=  1 
2 


(41) 


Then  we  have 


(*+-) 

a,c^H{zW  2  ) 

2  2 


0^it< 


M 


—  ^k<M 
2 


F^'\z)  =  a]c^H{zW^'*^^)  + 


0<it< 


M 


(42) 


and 
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pW(r)  =  ar<*^(r)/r<*>(0).  (45) 

Fig.  6  shows  the  signal  correlation  coefficient  for  the  different  subbands.  The  FIR  coefficients  are 
computed  based  on  [21],  and  the  number  of  taps  is  48,  64,  128,  and  256  for  A/=  4,  8,  16,  and  32, 

respectively.  From  the  figure  it  is  clear  that  is  flat  over  the  subbands  except  the  minor 

M  M 

degradation  at  the  two  sub-subbands  of  A:= - 1  and  —  when  r/Jy  is  not  an  integer.  Therefore, 

*  2  2 

the  signal  correlation  is  almost  equally  enhanced  for  the  different  subbands  after  the  subband  signal 
processing  based  on  modified-QMF  filter  banks. 

The  signal  correlation  coefficient  is  plotted  as  the  dashed  line  in  Fig.  3  for  Bx  -  1  and  k=0.  It  is 
evident  that  the  modified  QMF-based  filter  banks  have  very  close  signal  correlation  enhancement 
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performance  to  that  of  the  ideal  FIR  filters. 

In  the  above  analysis,  M  depends  on  the  signal  correlation  improvement  and  the  requirement  of 
subband  processing.  In  practice,  M  should  be  chosen  depending  on  the  environment  coherence 
bandwidth  [22]  and  the  employed  signal  bandwidth. 

IV.  Subband  Adaptive  Array 

A.  Concept 

In  the  above  sections,  we  have  discussed  the  signal  equalization  at  a  frequency-selective  fading 
environment  as  the  signal  correlation  between  multipath  rays  is  increased,  and  the  signal  correlation 
is  enhanced  by  subband  signal  processing  via  exploitation  of  the  filter  banks.  Therefore,  applying  the 
subband  signal  processing  technique  to  adaptive  arrays  forms  subband  adaptive  arrays,  which  is 
possible  to  solve  both  of  the  CCI  and  ISI  problems  in  a  multipath  fading  environment.  In  this  section 
we  discuss  the  concept  and  the  behavior  of  subband  adaptive  array. 

We  now  consider  the  adaptive  array  implementing  the  Wiener-Hopf  solution,  such  as  the  LMS 
adaptive  array,  with  the  modified  pseudo-QMF  filter  banks  here.  The  output  of  each  array  sensor  and 
the  reference  signal  are  divided  into  M  subbands,  respectively.  The  array  is  near-perfectly 
reconstructed  after  converging  to  the  reference  signal.  The  steady  state  optimum  weight  vector  at  the 
Arth  subband  is  thus  given  by 

k=\,...,M  (46) 

where 

=£[x<*>(0x‘*>''(/)]  (47) 

ri,*>=£[x<*>(0r<*’'(0]-  (48) 

From  the  discussion  in  Sections  2  and  3,  the  multipath  rays  in  each  subband  encounter  increases 
signal  correlation,  and,  asymptotically  they  behave  like  a  single  ray  for  high  values  of  M 
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B.  Performance  analysis 


From  (8)  and  (12),  the  array  output  signal-to-inierference-plus-noise  ratio  (SINR)  without  subband 
signal  processing  is  obtained  as, 


S1NR„  = 


Pso  ^  } 

“  P,,,  £{w''R^.w} 

-1—  “  ^£>^Z) 


(49) 


where  ^so  is  the  SOI  output  power  and  the  output  interfering  signal  power,  respectively.  With 
subband  signal  processing,  the  SOI  output  power,  ,  and  the  output  interfering  signal  power, 
PjP  ,  at  the  Ath  subband  are  given  by 


(50) 


and 


pW  _ 


~°D  ~S  “D 

respectively,  where 

g<*>  =  cr<*>^(l-crra'f>''RW-'aif>). 

The  output  SINR  of  a  subband  adaptive  array  is  expressed  as 


(51) 


(52) 


SINR„„=-i=^ 


Ypp  2:g<*>vrarRrsif> 


SUB  A/-1 

Y.C' 

k=0 


^'^(*)23(A)^rOhjo) 

k^O 


(53) 


The  SINR  improvement  by  subband  array  processing  is  most  significant  in  the  following  cases. 
The  first  case  is  that,  the  number  of  DOF’s  is  less  than  the  total  number  of  rays,  but  is  larger  than  the 
total  number  of  the  desired  and  CCI  signals,  while  the  signal  correlation  between  the  multipath  rays 
is  small.  Consider  a  mobile  station  at  high-speed  digital  communications  environment.  Although  the 
number  of  CCI  signals  is  limited,  they  can  arrive  with  a  large  number  of  multipath  rays  with  low 
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signal  correlation.  In  this  case,  the  conventional  adaptive  array  cannot  sufficiently  suppress  the  CCI 
rays  as  well  as  the  interfering  multipath  rays  of  the  desired  signal,  whereas  the  subband  adaptive 
array  can  suppress  the  CCI  rays  and  combines  the  multipath  rays  of  the  desired  signal. 

Another  case  is  that,  the  multipath  rays  of  the  desired  signal  come  from  a  small  angular  sector  so 
that  it  is  not  possible  for  conventional  adaptive  array  to  suppress  the  interfering  multipath  while 
maintaining  the  proper  gain  of  the  desired  signal.  A  subband  adaptive  array  is  not  sensitive  to  the 
problem  since  it  is  not  necessary  for  the  array  to  suppress  the  multipath  rays. 

V.  Performance  Simulation 

Computational  simulations  are  performed  to  confirm  our  analysis.  An  equi-spaced  linear  3- 
element  adaptive  array  is  considered  and  the  interelement  spacing  is  half  wavelength.  Sample  matrix 
inversion  (SMI)  method  [23]  is  used  to  approximate  the  optimum  weights  in  obtaining  the  output 
SINR  and  bit  error  rate  (BER).  Both  the  desired  signal  and  the  interference  signals  are  considered 
with  QPSK  modulation,  and  the  roll-off  rate  is  0.5. 

We  first  consider  the  case  where  only  the  desired  signal  with  two  rays  is  present.  It  is.  assumed  that 

ah,  =  cxhi  -  dB,  and  cr^  =0  dB.  Fig.  7  shows  the  output  SINR  versus  the  number  of  subbands, 
M  (In  the  figures  hereafter  M=1  denotes  the  case  where  the  signal  is  processed  without  subband 
decomposition).  In  Fig.  7,  the  DOA  of  the  direct  ray  is  0  degree  from  the  broadside  direction.  We 
show  two  cases  where  the  DOA  of  the  delayed  ray  is  from  40  degrees  and  2  degrees,  respectively. 
The  time  delay  between  the  two  rays  is  so  that  the  two  rays  are  uncorrelated. 

For  the  case  of  =40  degrees,  the  two  rays  are  widely  separated  in  space.  With  two  DOF’s,  the 
uncorrelated  delayed  ray  is  suppressed,  leading  to  high  output  SINR  is  high,  even  when  no  subband 
signal  processing  is  performed.  Thus,  the  improvement  by  subband  signal  processing  is  not 
significant.  However,  for  the  case  of  ^c:=2  degrees,  the  output  SINR  is  very  low  (1.7dB)  when  no 
subband  signal  processing  is  performed.  The  degradation  of  the  output  SINR  is  caused  by  the 
resolution  limitation  of  the  array  since  the  two  signals  are  closed  separated.  The  array  output  SINR 
performance  is  improved  rapidly  as  M  increases,  which  confirms  the  equalization  ability  of  the 
subband  signal  processing. 
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We  add  other  two  multipath  rays  into  the  above  case.  The  corresponding  array  output  SINR  is 
shown  in  Fig.  8.  The  new  multipath  rays  have  the  same  input  signal  level,  and  the  DOA's  are  -20 

and  -60  degrees,  and  their  delay  times  relative  to  the  direct  ray  are  2  7^  and  3  ,  respectively.  In 

this  case,  when  no  subband  signal  processing  is  performed,  the  four  rays  are  uncorrelated  to  each 
other,  and  the  array  will  try  to  suppress  the  three  delayed  rays.  We  note  that  the  3-element  array  only 
has  two  DOF’s.  Therefore,  the  array  DOF’s  are  not  sufficient,  and  the  array  cannot  provide  good 
performance.  It  is  clear  from  Fig.  8  that  the  output  SINR  is  near  OdB.  On  the  other  hand,  when 
subband  signal  processing  is  performed,  the  four  rays  will  asymptotically  behave  as  an  equivalent  ray 
and  the  output  SINR  increases  rapidly  as  M  increases.  The  output  SINR  exceeds  13dB  when  A/ =  4 

or  larger,  and  the  BER  is  lower  than  1 0"^ . 

Fig.9  depicts  the  output  SINR  when  two  CCI  signals  are  added  to  the  case  considered  in  Fig.7. 
The  CCI  powers  are  =  aj,  =20dB,  and  their  DOA’s  are  -20  and  -60  degrees,  with  delay  time 
difference  equals  to  .  When  no  subband  signal  processing  is  performed,  the  array  will  suppress  the 
delayed  ray  of  the  desired  signal  and  the  two  uncorrelated  CCI  rays.  Similar  to  the  case  in  Fig.8,  the 
number  of  DOF’s  is  not  sufficient  to  achieve  this  task.  The  result  is  very  poor  output  SINR  as  worse 
as  -1.6dB.  When  subband  signal  processing  is  performed,  the  two  rays  of  the  desired  signal  will  be 
treated  as  an  equivalent  single  one,  and  the  two  rays  of  the  CCI  signal  as  another  equivalent  ray. 
Thus,  the  array  DOF’s  will  be  sufficient,  and  the  output  SINR  increases  to  7.1dB  for  M  =  4,  and 
12dB  forM=  8.  The  corresponding  BER  performance  is  shown  in  Fig.lO,  which  decreases  rapidly  as 
M increases. 


VI.  Conclusion 

We  have  pr6posed  subband  adaptive  arrays  for  multipath  fading  mitigation  in  wireless 
communications.  It  has  been  shown  that  subband  adaptive  arrays,  which  provide  sub-optimal 
frequency-spatial  signal  processing,  are  very  effective  both  in  mitigating  the  inter-symbol 
interference  (ISl)  problem  caused  by  frequency-selective  fading,  and  in  suppressing  the  co-channel 
interference  (CCI)  signals  in  mobile  communications. 

The  performance  of  the  adaptive  array  in  the  presence  of  multipath  rays  has  been  analyzed.  We 
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have  shown  that  these  rays  can  be  considered  as  a  single  path  when  the  signal  correlation  between 
the  rays  is  high.  The  analysis  of  the  signal  correlation  enhancement  between  the  multipath  rays  by 
subband  signal  processing  was  presented.  We  have  discussed  the  signal  correlation  enhancement 
performance  for  filter  banks  using  ideal  filters,  DFT,  and  modified  QMF  filters.  The  excellent 
performance  of  modified  QMF  filter  banks,  which  are  practical  to  use,  is  also  demonstrated.  The 
simulation  results  confirmed  the  analyses  and  the  effectiveness  of  the  proposed  method  to  combat  the 
multipath  fading  for  both  the  desired  signal  and  the  CCI  signals. 
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Fig.  1  Block  diagram  of  a  filter  bank 
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Fig.4  Signal  correlation  coefficient  of  DFT  filter  banks  at  the  different  subbands 
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Fig.7  Output  SINR  vs.  M 
(<Toi  =<702=10^^®*  Qdi  =0deg.,  T=r5) 
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Fig.8  Output  SINR  vs.  M 

(^DI  =  0-03  =  lOdB.  e^,^  =Odeg.,  6^,^  =20  deg., 
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Fig.9  Output  SINR  vs.  M 
(CT^ji  ~  O’ 1)2  —  lOdB,  0^1  —Odeg.,  0p2  “* 
(^/i  =^/2  =20dB,  =-20  deg.,  Oj^  =* 


20  deg.,  T^=T5, 
-60  deg.,  rj^Ts) 
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